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Abstract 

> , 

We solve a long-standing problem in particle physics: that of deriving the 
Deep Inelastic structure functions of the proton from the fundamental theory 
of strong interactions, Quantum ChromoDynamics (QCD). In the Bjorken 
limit, the momenta of the constituents of the proton (the partons) can be 
assumed to be in a two-dimensional plane in Minkowski space: a dimensional 
reduction of QCD to two space-time dimensions. Two dimensional QCD is 
then shown to be equivalent for all energies and values of number of colors N 
to a new theory of hadrons, Quantum HadronDynamics (QHD). The phase 
■ space of QHD is the Grassmannian (set of subspcaes) of the complex Hilbert 

space L 2 (R) . The natural symplectic form along with a hamiltonian define a 
classical dynamical system, which is equivalent to the large N limit of QCD. 
't Hooft's planar limit is the linear approximation to our theory: we recover 
his integral equation for the meson spectrum but also all the interactions of 
the mesons. The Grassmannian is a union of connected components labelled 
by an integer (the renormalized dimension of the subspace) which has the 
physical meaning of baryon number. The proton is the topological soliton: 
the minimum of the energy in the sector with baryon number one gives the 
structure functions of the proton. We solve the resulting integral equations 
numerically; the agreement with experimental data is quite good for values 
of the Bjorken variable xb > 0.2 . 
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Chapter 1 
Introduction 



The first indication that the atom contains a point-like nucleus came from 
experiments of Rutherford. He scattered alpha particles (which are positively 
charged) from a gold foil. Since the negative charges (electrons) were known 
to be much lighter than the alpha particles, they would not scatter the alpha 
particles very much. If the positive charges inside the atom were more or less 
uniformly distributed, the alpha particles would not be scattered through 
wide angles. Rutherford found to the contrary that the alpha particles were 
scattered by wide angles. In fact the probability of scattering through an 
angle 9 is proportional to sin -4 | , exactly what would happen if all the 
positive charge of the atom were concentrated at a point: he had discovered 
the atomic nucleus. Soon after, it was realized that the nucleus is of finite 
size, although small compared to the atom. It is composed of protons and 
neutrons bound together by the strong interaction. 

Many years later, another celebrated series of experiments studied the 
structure of the proton itself. It was found that it was not an elementary 
particle either, in fact that it was made of point-like constituents. In such 
a 'Deep Inelastic Scattering' experiment an electron (or neutrino) beam is 
scattered by a proton (or a nucleus containing both protons and neutrons). 
The idea is to learn about the still mysterious strong interactions using the 
electroweak interactions as a probe. 

A particle that can take part in the strong interactions is called a 'hadron'. 
There are two kinds of hadrons: those of half-integer angular momentum are 
called 'baryons' and those of integer angular momentum are the 'mesons'. 
The proton is the lightest baryon and the 7r -meson (pion) the lightest 
meson. There are an infinite number of baryons and mesons but the more 
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energetic ones are unstable against decay into the lighter hadrons. (We con- 
sider an idealized world in which the electromagnetic, weak and gravitational 
interactions of the hadrons with each other are ignored. Thus a particle is 
considered stable if its decays are purely electromagnetic or weak. This is 
the sense in which the pion is stable. ) 

The hadrons are composed of more elementary constituents: the quarks, 
anti-quarks and gluons. These constituents are collectively called partons 

However, it has been found experimentaly that it is impossible to create 
the partons in isolation: they only exist inside hadrons, and thus, are not 
true particles. This is the phenomenon of 'confinement'. More precisely, 
only combinations of quarks and gluons invariant under the action of the 
color group SU(N) (see below) exist as isolation. Thus hadrons can be 
defined to be states that are invariant under color. Some authors use the word 
hadron to include the quarks and gluons. We will use the term 'hadron' to 
refer to a bound state of quarks and gluons which can exist as an isolated 
particle. Of course such a particle may be unstable against decay into other 
hadrons. 

Around the same time of these developments an entirely different picture 
of a baryon was proposed by Skyrme ||: that it is a topological soliton 
made of an infinite number of mesons. Unfortunately, this idea did not fit 
with the dominant view of the time and was largely ignored. In the mid- 
eighties Skyrme's idea was finally revived by a group at Syracuse University 
(including the author) J|, |, || and integrated into the modern theory of 
strong interactions, i.e., Quantum Chromodynamics. Witten |7] clarified 
why the baryons are fermions (when N is odd) even though the underlying 
theory has only bosonic fields. ( Witten had already arrived at the idea that 
baryons are solitons in the large N limit of 't Hooft || by independent 
arguments ||. ) 

But then the problem remains how to reconcile it with the picture of a 
baryon as a bound state of point-like constituents. We will show in these 
lectures how the parton model can be derived from the soliton model. It fact 
it had never been possible to derive the distribution functions of the partons 
inside the hadron from fundamental principles. Our picture will solve this 
problem, providing for the first time a quantitative theory of the structure 
of a proton. 

The basic idea of our approach is to find a new description of strong in- 
teractions directly in terms of hadrons rather than in terms of quarks and 
gluons. It should be equivalent to the color singlet sector of QCD at all ener- 
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gies and all values of the number of colors. We will call this theory quantum 
hadrondynamics. In four dimensions this is still just an idea; we havent been 
able to show yet that such a new paradigm for strong interactions exists. 
However in the two dimensional case, I constructed such an equivalent alter- 
native formalism some years ago. Part of the motivation for studying the two 
dimensional case was that it would provide a solution to the problem of de- 
riving the structure functions of four dimensional QCD. (It has been known 
since the early work of Feynman and Bjorken that Deep Inelastic Scattering 
can be explained by a two dimensional theory of strong interactions.) The 
first report of the idea of such an equivalent theory was in Ref. The 
solutions of the integral equations for the baryon wavefunction was first stud- 



ied in Ref. ||1 1|| . A more detailed description of the theory appeared in Ref. 



12]. In another direction, these ideas were applied to spherically symmetric 



situations in Ref. |13 . The methods of geometrical quantization were ap- 



plied to the problem in Ref. [H| . There a version with bosonic quarks was 



also studied. We returned to the study of the baryon wavefunctions in Ref. 
T5| where the derivation of the parton model, already mentioned in JT2] 



was 



given in more detail. In papers which will soon appear JT6[ we will extend 
these results to include the anti-quark and sea quark distributions. Some of 
the ideas that go into this work are outlined in the last chapter. The basic 
mathematical ideas necessary to derive the gluon structure functions have 
also been developed [|17]]. We hope to solve that problem as well in the near 



future. 

The collinear approximation to QCD has alos been studied by other pa- 
pers [0 However the combination of these ideas with the picture of a baryon 
as a soliton as well as the idea of a quantum hadrondynamics seem unique 
to our approach. 

The ideas we introduce range from infinite dimensional geometry to the 
phenemenology of particle physics. Indeed a whole new set of tools had to 
be developed in order to implement the new paradigm we propose. Some 
of the methods are currently not associated with particle physics: they are 
more closely related to classical mechanics and are perhaps more familiar to 
mathematicians. But the new ideas of one generation become the standard 
lore of the next. I hope these lecture notes will prepare a new generation of 
theoretical physicists to pursue these ideas. Prior knowledge of quantum field 
theory is not essential, but will provide perspective. In general, a higher de- 
gree of mathematical maturity than specialized knowledge of particle physics 
is assumed. 
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In the next few sections of this introduction we give a quick summary of 
the physics background necessary. The reader who finds this boring should 
skip to the next chapter, after a glance at sections 1.4 and 1.5 which are 
essential to the argument. 



1.1 Deep Inelastic Scattering 

Let us consider in more detail the scattering of an electron by a proton. The 
electron emits a virtual photon of momentum q = kf — ki where ki and 
kf are the initial and final electron momenta. This photon then interacts 
with the hadron producing some state \X > which could contain many 
particles: several mesons and some excited baryons and anti-baryons. This 
is an inelastic scattering process. If we sum over all the states \X > so 
produced, we will get the 'inclusive' cross-section for the scattering of the 
electron. Upto well-known electromagnetic effects, the cross-section for this 
scattering is given by the matrix element [ITS 



W, u (P,q) = ^Y,Y,< P °W)\ X >< X \M0)\P°>(^) 4 5\P + q-P x ) 



<7r - x 



l -Y^J dxe-*** < Pa\J,(x)J u (Q)\Po > 



Here, Jfi(x) is the electromagnetic current operator which is the one 
appropriate for electron scattering: the weak interactions of the electron are 
much smaller in comparison. For neutrino scattering it would instead be the 
charged weak current (correponding to W ) if the neutrino is converted into 
an electron; and the neutral current (corresponding to Z° ) otherwise. Also, 
P is the momentum of the hadron and we assume that we have averaged 
over the spin o of the hadron. (The target is usually not polarized, so 
that we must average over the values of the spin variables. More detailed 
measurements with polarized targets have been made more recently; we will 
not consider these for now.) 

The tensor is symmetric and transverse: conservation of the electric 

current gives 

q»W, u (P,q)=0 (1.1) 
Lorentz invariance and parity (which is a symmetry of the electromagnetic 
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interactions) imply that the tensor W^ip, q) has the form 



W,„(P,q) 



W 1 (x B ,<?) + (1.2) 

Being Lorentz scalars, the 'structure functions' Wi j2 can depend only 
on the Lorentz scalars q 2 and P ■ q ( P 2 is fixed to have the value 
m 2 , where m is the mass of the proton). When an electron is scattered 
against a target, the momentum q of the photon it emits is space-like. It is 
conventional to use the positive number Q 2 = —q 2 to describe the energy 
of the photon. The higher the value of Q 2 the better the resolution of 
our measurement of the structure of the hadron. As the other independent 
variable we can take the dimensionless ratio ('Bjorken variable') 

If the hadron we are studying is stable against strong decays (which is always 
the case for experimental reasons ) it will be the lightest particle with its 
quantum numbers. Then, the mass of the intermediate state will be greater 
than or equal to the mass of the target: 

(P + qf >m 2 ; [ 1] > 0. (1.4) 

Xb 

This shows that the Bjorken variable takes values in the range 

1 > x B > 0. (1.5) 

It is often more convenient to use another equivalent pair of structure 
functions 

F 1 (x B ,Q 2 ) = W 1 (x,Q 2 ), F 2 (x B ,Q 2 ) = ^-lw 2 (x B ,Q 2 ). (1.6) 

In its rest frame the hadron has a certain 'size' a : it is of the order of 
the the charge radius of the proton, a -1 ~ lOOMeV . If we take the limit 
Q 2 » a~ 2 keeping x B fixed we are looking deep inside the hadron: this is 
the region of Deep Inelastic Scattering. The basic idea is much like that of a 
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microscope: to see inside an object of size a we need light of a wavelength 
that is small compared to a ; or equivalently an energy for the photon that 
is large compared to a -1 . 

When the states are normalized by the usual convention 

< p\p' >= 2E p (2tt) 3 5 3 (P - P') } (1.7) 

the functions Fi 2 {xb, Q 2 ) are dimensionless. It is found experimentally 
that these functions are approximately independent of Q 2 : they depend 
essentially only on the dimensionless ratio x b ■ The simplest explanation of 
this phenomenon is that the proton is made of pointlike massless constituents: 
the structure function would then depend only on the dimensionless ratio xb 

If we multiply P ■ q and Q 2 by the same number, the structure 
functions are approximately invariant: this is the approximate symmetry 
of 'scale invariance'. In the first versions of this model, the partons were 
concieved of as free particles. Of course they must interact to bind into a 
hadron, but the idea was that at high enough Q 2 the interaction would 
for some reason be small. This was the germ of the idea of 'asymptotic 
freedom': that the constituents of the hadrons behaved like free particles at 
short distances, or high energies. 



This 'parton' model was subsequently derived |2(J, p^j , p3|| from a 
much deeper fundamental theory of strong interactions, Quantum Chromo- 
dynamics (QCD). This is a non-abelian gauge theory which has the un- 
usual property that at high Q 2 the coupling constant (which measures the 
strength of the interaction) vanishes like — W . Here A is a parame- 

ter with the dimensions of momentum which is a fundamental parameter of 
QCD. Hence at high energies (compared to A ) QCD tends to a free theory, 
yielding asymptotic freedom and the parton model. In the limit of large Q 2 
, the structure functions are predicted to be independent of Q 2 . At finite 
values of Q 2 , the structure functions have a slow dependence on Q 2 ; as a 
polynomial in the coupling constant, or equivalently, — W . Moreover, we 

can calculate this dependence on — W using perturbation theory: scale 

invariance is broken in a way that is calculable. Perturbation theory allows 
us to calculate the Q 2 dependence of the structure functions: given the 
value of Fi^xb, Q 2 ) at one (large) value of Q 2 , we can calculate it for 
any other large value of Q 2 . ( We mean that Q should be large compared 
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to A , which is of the order of 100 MeV .If ^ is not large, perturbatiion 
theory can no longer be used.) 

It is one of the triumphs of modern particle physics that perturbative 
QCD accurately describes these scale violations. 

However, the Xb dependence of the structure functions are not as well- 
understood. Perturbation theory is not sufficient to calculate it: F 2 (xb, Q 2 ) 
is in a certain sense the probability distribution (in momentum space) of a 
parton inside a hadron. Thus it describes how partons bind together to form 
a hadron. Perturbation theory around a free field theory can never describe 
such bound states. Due to the lack of any fundamental understanding of the 
xb dependence of the structure functions, physicists have been forced to ex- 
tract them directly from experiment. A whole generation of experimentalists 
and phenemenologists have worked to produce a quite reliable extraction of 
the xb dependence of structure functions from data. It will be our goal to 
explain this from the fundamental theory of strong interactions, QCD. 

In the next section we will summarize the definition of QCD as a quantum 
field theory. 



1.2 Quantum Chromodynamics 

Quantum Chromodynamics is the fundamental theory of strong interactions. 
It is a non-abelian gauge theory with gauge group SU (N) with matter fields 
(quarks ) which are spin half fermions transforming under ( Nf copies of 
) the fundamental representation of SU(N) . The natural number A" 
is called the 'number of colors'. It has the value 3 in nature; but it 
will be convenient to leave it unspecified until the very end when we make 
comparisons with data. The 'number of flavors' Nf is 2 for most purposes 
although it can be as high as 6 in principle: the heavier flavors of quarks 
can be ignored for most purposes. Again, it is best to leave Nf as an 
arbitrary parameter for now. 

The action principle that defines QCD is 

N r Nf f 

S=— J trF^F^d 4 x + J2J <f H7 • V + m a ]q a d 4 x (1.8) 

(2=1 

The Yang-Mills field strength is 

Ffj, u = d^A v - d v A^ + [Ap, A u ] (1.9) 
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where are a set of four N x N anti-hermitean matrices: A^dx 11 

is a one-form taking values in the Lie algebra of the unitary group U(N) 
. In perturbation theory, this 'gauge field ' describes 'gluons': massless spin 
one particles analogous to the photon. q a is a fermionic field (describing 
quarks) which is a Dirac spinor transforming in the fundamental represen- 
tation of the 'color' group U(N) , for each a . The index takes values 
a = 1 • ■ ■ Nf , where Nf is the number of copies of such quarks: the number 
of flavors. 

The parameters of the theory are (in addition to N and Nf ) , the quark 
masses m a for a — 1, • • • Nf and the dimensionless 'coupling constant' 
a . In fact the quantum theory is defined by an additional set of rules for 
'renormalization' which replace the constant a by a new parameter A 
which has the dimensions of mass as well. These issues are well-known and 



have been reviewed in several articles 23 , so we will avoid them here. T he 



basic point is that the the coupling constant acquires a dependence on the 
energy scale a ~ — W . 

These Nf + 1 parameters should in principle determine the masses and 
decay rates as well as the structure functions of all the hadrons. It has turned 
out to be quite difficult to predict hadronic properties from this fundamental 
theory. We need to develop approximation methods which make the problem 
tractable. For the most part this problem is still not solved. In these talks 
we will describe how we solved a part of this puzzle: that of calculating the 
hadronic structure functions from QCD. 



1.3 Structure Functions from QCD 

We noted earlier that since the electromagnetic interactions are well-understood, 
it is sufficient to concentrate on the 'unknown' part of the problem, which 
is encoded in the structure functions Wi^xbiQ 2 ) ■ These were defined in 
terms of the expectation values of the product of current operators. Now, 
part of QCD is tractable by perturbation theory; this allows us to reduce the 
'unknown' part of the problem further to the expectation values of simpler 
operators: the structure functions can be expressed in terms of 'parton dis- 
tribution functions'. The problem we will solve is that of calculating these 
distribution functions. 

The basic result that makes this reduction possible is called the 'factor- 
ization theorem': the structure function can be written as a convolution of 
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two factors: one whose xb and Q 2 dependence can be calculated in 
perturbation theory (the 'hard' factor); and another (the 'soft' factor) which 
has an unknown xb dependence but which is independent of Q 2 . This is 
a result of quantum field theory that can be proved to all orders of perturba- 
tion theory; so it is comparable in depth to the proof of renormalizability of 
gauge theories. We will just give an intuitive argument, pointing the reader 



to the literature on perturbative QCD for details |L9|, £4 . 

The photon is scattered by a quark of charge e a inside the hadron. To 
leading order, the probability of the photon being scattered by the proton 
can be written as the product of two pieces: the probablity that a quark 
of some momentum k will scatter the photon and the probability that 
there is such a quark inside the proton. (Then we of course sum over the 
intermediate momentum k .) This is the simplest version of the factorization 
theorem. In the next order of perturbation theory, we will have to include the 
possibility that quark might emit a gluon before absorbing the photon and 
then re-emit it. We can imagine the scattering of the photon by the quark as 
a subprocess with its own structure function, except that it can be calculated 
within perturbation theory. The probability that a hadron contains a quark 
of a given momentum is of course not calculable in perturbation theory: that 
is the distribution function which we will attempt to understand. 

In other words, the structure function of the proton can be written as 
the convolution of two pieces: the structure function of a parton (which is 
computable perturbatively) and the distribution function of the parton inside 
the hadron (which is non-perturbative). 

In more detail, 




(1.10) 



Here <fi a {%B) is the probability of finding a quark of flavor and momentum 
fraction xb inside a proton; ip a is the probability for an anti-quark of 
flavor a , and (pc f° r gluons. Note that all the Q 2 dependence is 
in the first factor; it can be calculated as a power series in — W using 

l°g 7^ 

the standard rules of perturbation theory. To leading order in perturbation 
theory, even the first factor is independent of Q 2 : 

C«(x B , Q 2 ) = Cl(x B , Q 2 ) = e 2 a 8(x B - 1), C?(x B , Q 2 ) = 0. (1.11) 
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If bi(k), b\(k) are the creation-annihilation operators for the parton, we 
have 

U*b) = E/ |^ < P(j\b\{x B PM)bi{x B PM)\P(T > . (1.12) 

Here, we are averaging over all possible values of the part of the momentum 
orthogonal to the plane spanned by P and q , (called the 'transverse 
momentum ). 

In terms of the field operators, we get the quark distribution functions 

<t>a{x B ) = ^E / dye~ lXBPy < Pa\q a {y^ T )l-qamP<J > (1.13) 

^ a J 

and the anti-quark distribution functions 

M*b) = —£,1 dye' lXBPy < Pcr\q a a (y,O T )b-haq a p(0)\PcT > ■ (1.14) 

The position arguments of the field operators are separated by a null line in 
the plane spanned by P and q . The averaging over transverse momenta 
implies that in position space, these field operators have the same transverse 
co-ordinates. 

Since the gluons do not carry electric charge, to the leading order of 
perturbative QCD, the gluon distribution functions are not necessary in order 
to understand the structure of the proton. We will ignore them for now. 

These formulae are not gauge-invariant and are to be understood in the 
null gauge A_ = . In a general gauge we should insert a parallel transport 
operator PeJ A -( y, ° T ^dy along the null line connecting two field operators. 



1.4 Reduction of QCD to Two Dimensions 

We are probing the hadron (whose momentum is a time-like vector P ) 
by a photon (or W, Z -boson) with a space-like momentum q . These 
two vectors define a two-dimensional subspace of Minkowski space. There 
is also a space-like vector a which characterizes the 'size' of the hadron: 
the length of a is the charge radius of the hadron. In the Deep Inelastic 
limit, Q 2 = —q 2 » | a, | 1 keeping xb = ^p— fixed. This means 
that the hadron is of very large size in the directions transverse to the two 
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dimensional subspace spanned by P and q . By the uncertainty principle, 
the transverse momentum of its constituents is of order | cz | 1 , which is thus 
small compared to Q . A reasonable first approximation would be to let 
the hadron have infinite extent in the transverse directions; and to require 
the momenta of the constituents to lie entirely in the plane spanned by P, q 
. In other words the approximation is that the fields are independent of the 
transverse spatial co-ordinates. 

This is analogous to the procedure of dimensional reduction popular in 
unified field theories of gravity. The main difference is that the dimensions 
that are ignored do not form a subpace of small volume: instead they are 
infinite in extent. The point is that the momenta in these directions are 
small which is the same as requiring that the fields are independent of those 
spatial directions: a dimensional reduction. 

Thus Deep Inelastic Scattering is described by the dimensional reduction 
of QCD to two dimensions. That is what makes this phenomenon accessi- 
ble: two dimensional gauge theories have been understood by the large N 
method. With some further (less drastic) approximations we will be able to 
determine the spectrum and structure of hadrons in two dimensions. These 
should then be verifiable experimentally. 

The dimensional reduction of QCD to two dimensions is given by the 
action principle 

N r Nf r 

S = — trF^F^x + ]T / ^[-i-y • V + m a ]q aa d 4 x 
4«i J a=1 J 

N f N f f 

+ 2^J t ^,<P A ) 2 d 2 x + —J tr[0 3 ,0 4 ]Va; + y q aa {-i)T A a ^ A q aP d 2 x. 

Here we allow the indices /i, v to take only the values 0, 1 . The gauge 
field now splits into a 1-form in the two dimensional space with components 
and a pair of scalar fields 03 = A 3 , 4 = A 4 corresponding to the 
transverse polarization states of the gluon. 

The four dimensional Dirac spinor q splits into a pair of two-dimensional 
spinors q a , corresponding the two eigenvalues of 7374 . Moreover a\ is 
the coupling constant of the two-dimensional theory which has dimensions of 
(mass) 2 . It is a combination of a and the size of the hadron: a x ~ cra~ 2 
. Thus the reduction of QCD to two dimensions gives a gauge theory with 
scalar fields in the adjoint representation, and twice as many flavors of quarks 
as the original theory. 
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We will study for the most part two dimensional QCD, defined by the 
action 

N r Nf r 

S = — / trF^F^d 2 x + J2 q aa Hl-V + m a ]q aa d A x. 

This is not quite the same as the dimensional reduction of four dimensional 
QCD to two dimensions: there are no scalar fields. This truncated theory will 
be sufficient to to determine the quark and anti-quark structure functions to 
the accuracy we need. The essential techniques required to solve the theory 
including the scalar fields have been developed as well. We will return to 
them in a later publication. 



1.5 Two Dimensional QCD in Null Gauge 

The essential simplification of two dimensional QCD is that the gauge fields 
can be removed completely from the problem, leaving just the quark fields 
q and q and the scalar gluons <pA ■ This should not be too surprising: in 
D dimensional space-time, a gauge field has D — 2 polarization states. ( 
Of course, we still have the scalar fields which are the remnants in two 
dimensions of the two polarization states of the gluon in four dimensional 
gauge theory. But we will ignore them, as they are mostly relevant to the 
determination of the gluon structure functions: a problem we will postpone 
to a later publication.) 

We will change notation slightly from the last sectiomthe indices a, b 
etc. will denote the pairs aa , b(3 , so they will range over 1, 2, • • - 2Nf 
.This simply reflects the fact that the two dimensional theory has twice as 
many flavors as the four dimensional theory: the transverse polarization label 
looks just just like a flavor index to the two-dimensional theory. 

In the null gauge, the action of the theory becomes (see Appendix A) 

{ 1 777 2 1 N 

-\p+ — } -iA t jxa + ^ tr[d x A t ] 2 . (1.15) 

The field A t carries no dynamical degrees of freedom: it can be elimi- 
nated in terms of x by solving its equations of motion: 

= / \\ x - v\ ■ x ]ai {y)xa 3 {v) ■■ dy. (1.16) 
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(The Coulomb potential in one- dimensional space is given by ^\\x — y\ 
S(x — y) .) The hamiltonian of the resulting theory is 



H = fdx X ^-[p+— }Xai 
J 2 p 

1 r 1 

iV^ 1 J 2^ : x]m ^ Xa ^ x "> X^(y)Xbi(y) ■ dxdy 



2N 

Now define the operator 

K(x,y) = -^:x Wt (x)xm(y): (1.17) 

which is gauge invariant. It describes the creation of a quark at x and 
an anti-quark at y , but in a color invariant combination: in other words 
it describes a meson. The hamiltonian above can be expressed entirely in 
terms of this operator after some reordering of factors of x ■ 

We can rearrange the quartic operator in the x ' s as a quadratic operator 
in M . First, 

: X iai (x) Xaj (x):: X ^(y) Xbi (y): = : X^{x) Xaj {x) X ^ {y) Xhi {y) : 

+ :X 1ai (x)xdy):<0\x aj (x)x W (y)\0> 

+ < o\x Wt (x)xM(y)\o >: xa,^)x t6j (y) : 

Now, (see Appendix A) 

< 0\x 1m (x) Xbl (y)\0 >= N5 a b f ^ = N 8l l -[5{x -y) + e(x - y)\ 

where, 

e( x - y) = V ( sgn (p)e ip(x - y) ^-. 

J 2ir 

Thus, 

■■X ]ai {x) Xaj {x):: X ^{y) Xbi {y): = : X ^{x) Xaj {x) X ^ {y) Xbi {y) : 

+ : X tai (^)x«(y) :: y - y) + e(x - y)] 
+—[S(x -y) + e{x - y)] : Xai{x)x ]ai {y) ■ 
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On the other hand, 

2^ 2 



K(x,y)M b a (y,x) = (-) : x*«(x) X »(y) " \ '"i//) \„ ; (.r) : 



i 2 

ViVy 

The terms corresponding to other orderings of x will involve e(y, y) which 
should be interpreted as zero. So, 

/ N\ 2 - 

: \ •"(,•) \„ ; (,•!:: \ '•'( //) \ ,„(//) : = -(-J M« (x , y) M b a (y , x) 

+y[*(x-y) + c(x-y)](-y)<(y,x) 

Thus, 

N\ r 1 fi 2 ~ a dp 



i r i r /iv\ 2 - 

-«i y (^J MZ(x,y)M b a (y,x) 



2N 

N 2 



+ (- —^)K{x,y)e{x -y^dxdy 

We have dropped terms involving 5(x — y) \x — y\ since this product is just 
zero; moreover the two terms involving e(x — y) are equal and have been 
combined. 

Now we use the identity 

/\ \ r \ ~ dp 
-\x - y\e(x - y)f(x - y)dxdy = V / -f{p,p)—. (1.18) 
I 7T J V LlX 



Thus 



H 1 f 1 . /A ;~ a . . dp 



l ai J \\ x ~ y\^{x,y)M b {y,x)dxdy 



where 



7T 
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The commutation relations can also be expressed entirely in terms of the 
variables M : 

{M b (p,q),M d (r,s)} = 1 [s c b 2nS(q - r)[8 a d sgn (p - s) + M d (p, s)} 

-S a d 2n5(s - p) [51 sgn (r - q) + M d (r, g)]) . 

We notice that in the limit of large N , these commutators become small: 
it is a sort of classical limit. 

1.6 Constraint on the Variable M 

We have all the essential ingredients of a reformulation of two dimensional 
QCD in terms of the color singlet variable M , eliminating quarks and 
gluons from the picture. The commutation relations and the hamiltonian 
together inply the time evolution equations for M . However there is one 
more ingredient which we can miss at first: the set of allowed values of M 
. Being bilinear in the fermionic variables, M behaves much like bosonic 
variables: they satisfy commutation relations (rather than anti-commutation 
relations). Also, they have a classical limit; in our case this is the large N 
limit when there commutators become small. But their origin as fermion 
bilinears have a residual effect: they satisfy a quadratic constraint. This 
constraint is ultimately an expression of the Pauli exclusion principle; even 
in the large N limit the fact the underlying degrees of freedom are fermions 
cannot be ignored. 

Regard the classical variable p(x, y) obtained by taking the large N 
limit of jfX^( x )Xi(y) as a matrix in the variables x and y . It is in 
fact the density matrix of quarks. Its eigenvalues have to be between zero 
and one: zero when the state is completely empty and one when it is filled. 
For the color singlet states (hadrons), each state is either completely filled 
with quarks or completely empty. Hence the eigenvalue is either or one, 
and this density matrix is a projection operator: p 2 = p or, 

/ p( x , y)p(y, z)dz = p(x, z). (1.20) 

This is a quadratic constraint on the classical variable. 

Another (oversimplified) way to understand this constraint is that p(x, y) 
corresponds to a meson state where a quark is created at x and an anti- 
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quark at y ■ H we now create two mesons, with the position of the anti- 
quark of the first coinciding with the quark of the second, they will annihilate 
eachother leaving us with one meson!. 

We will usually use the normal ordered variable M(x, y) = 4 : X 1 ^ ( x )Xj(v) '■ 
, which has the advanatage that M = on the vacuum. In terms of it the 
constraint becomes, in matrix notation 

[e,M]+ + M 2 = 0. (1.21) 

M has to satisfy some technical conditions as well, but we will talk about 
them later. 

It is possible to prove these constraints by a straightforward calculation 
on the operators, even when N is held finite. ] (See the appendix of Ref. 
jl3| .) It is important to note that they only hold in the color singlet subspace 
of the fermionic operators. 

This means that it is possible to understand two dimensional QCD with- 
out ever mentioning quarks, antiquarks or gluons. The operator M is gauge 
invariant and can be the field variable of a theory that directly describes 
hadrons. From this point of of view quarks and gluons are just mathematical 
artifacts. This gives us a whole new paradigm for the theory of strong inter- 
actions. Instead of a 'Quantum ChromoDynamics" in terms of unobservable 
quarks and gluons but a new "Quantum HadronDynamics" in terms of the 
directly observable particles, the hadrons. 

In the next chapter we will develop just such a point of view. We will 
first present a classical theory which is equivalent to the large N limit of 
the above theory. Upon quantization we will recover two-dimensional QCD. 
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Chapter 2 



Two-Dimensional Quantum 
Hadron Dynamics 

In this chapter we will present first the classical hadron theory and then its 
quantization. We will see that the quantum hadron theory is equivalent to 
QCD in two dimensions. 

2.1 Grassmannians 

Let H = L 2 (R) (g> C N2f be the complex Hilbert space of complex- valued 
functions on the real lineQ. Define the subspace of functions whose Fourier 
transforms vanish for negative momenta: 

H+ = = jT^e**^} (2.1) 

and its orthogonal complement 

H- = = f_J^y PX ^\ ( 2 - 2 ) 

The functions in 7i + are boundary values of analytic functions on the upper 
half of the complex plane in the variable x . (Of course is related to 

the lower half of the complex plane in the same way.) Define the operator 

lr The number N%f is the number of flavors of the two dimensional theory; it is twice 
the number of flavors of the four-dimensional theory. 
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e ('the sign of the momentum') to be —1 on and +1 on H + 

Clearly 



e f = e, e 2 = 1. 



(2.3) 



The Grassmannian of H is the set of all its subspaces. To each subspace 
W , there is an operator $ which is —1 on W and +1 on W 1 - . 



Conversely any such operator corresponds to an orthogonal splitting of H 
: it will have eigenvalues ±1 and W can be identified as the subspace 
with negative eigenvalue. Thus the set of such $ 's is the same as the set 
of subspaces of Ji .The set of subspaces of 7i is the Grassmannian, which 
we can thus also regard as the set of operators $ satisfying the above 
condition. 

But we dont want to allow all such subspaces: since Ti is infinite 
dimensional, that would give a Grassmannian that is too big. For example, 
its tangent space will not be a Hilbert space, and it wouldnt admit a Riemann 
metric. We should only allow subspaces W that are 'at a finite distance' 
from . It is convenient to introduce the variable 



which measures the deviation from the standard point e , which corresponds 
to the subspace . We will require that M be Hilbert-Schmidt. That 
is, the sum of the absolute magnitudes squared of all its matrix elements is 
finite. (See Appendix B for a rapid summary of some functional analysis we 
will need.) This sum is a measure of the distance of W from H- . 
Thus we define the restricted Grassmannian to be, 



Gr(W, e) = {M f = M, [e, M]+ + M 2 = 0, trM 2 < oo}. (2.6) 



The set of all such M 's forms an infinite dimensional Hilbert manifold. 
This manifold is defined by a quadratic equation in the (real) Hilbert space 
of self-adjoint Hilbert-Schmidt operators. 



Again, 



(2.4) 



M = $-e 



(2.5) 



If we split the operator into 





(2.7) 
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we get (using the constraint on M ), 

a =-(btf + a 2 ), d= --(tfb + d 2 ), ab + bd = 0. (2.8) 

Since M as a whole is Hilbert-Schmidt, 6 : 7i + — > 7i_ is Hilbert-Schmidt 
as well. The above constraints then imply that a and d are trace-class. 
This means roughly speaking that the sum of their eigenvalues is absolutely 
convergent (see the appendix B for precise definition). 

Note that a is a positive operator and d a negative operator. In 
particular, tra > and trd < . 

The tangent space at the origin is given by the special case where a, b, c, d 
are all infinitesimally small. The constraints then show that a and d are 
second order infinitesimals while b are first order. In the first order the only 
constraint on b is that it has to be a Hilbert-Schmidt operator. Thus, the 
tangent space at the origin is a complex Hilbert space, consisting of operators 

of the form M — ^ ^ ^ ^ . In other words we have the identification, 

T Gr(H, e) = 1 2 {H + -> (2.9) 



2.2 The Dirac Theory of Fermions 

We have introduced the Grassmannian from a purely geometric point of view 
above. But it has a natural physical interpretation in terms of the Dirac 
theory of fermions. 

In Dirac's theory, the complex Hilbert space 7i represents the set of 
states of a single fermion. There is a self-adjoint operator h representing 
energy. In the null co-ordinates we are using, energy is 

Po=~\P+— ]■ (2-10) 
z p 

Thus it has the same sign as the null component of momentum p . The main 
physical obstacle to this interpretation is that the energy is not positive: all 
states are thus unstable with respect to decay into the negative energy states 
of lower and lower energy. 

Dirac's main idea was that in any physical state all except a finite number 
of the negative energy states are occupied. Since fermions obey the Pauli 
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exclusion principle, only one particle can occupy a given state: decay into 
states of very low energy is forbidden because those states are occupied. In 
particular the vacuum is not the state containing no particles. Instead the 
vacuum is the state where all the negative energy states are occupied and 
the positive energy states are empty. The filled negative energy states in 
the vacuum is called the 'Dirac Sea'. All physical quantities such as energy 
or charge are to be measured in terms of their departure from the vacuum 
value: even though the vacuum contains an infinite number of particles, it is 
still assiged zero energy, charge etc. 

An excitation from the vacuum could be a state where a finite number 
of positive energy states are occupied and a finite number of negative energy 
states are empty. Such a state will have positive energy compared to the 
vacuum. The empty negative energy states (also called 'holes') have positive 
energy. The holes themselves behave just like the particles except that they 
have the opposite value for some observables such as electric charge. They 
are the 'anti-particles'. The number of particles minus the anti-particles is a 
conserved quantity which can take any integer value. An arbitrary physical 
state is a linear combination of such states containing a finite number of 
particles and anti-particles. 

A mathematical interpretation of this situation can be given in terms of 
a modified exterior product of Ti . In familiar non-relativistic theories, the 
space of states of a multi-particle system of fermions would be the exterior 
power of order r , where is the number of fermions: J2^=o A r 7i . Instead 
in Dirac's theory we split the one particle Hilbert space into two subspaces 
H- © H + of negative and positive energy. The space of states of the multi- 
particle system is 

oo 

JF = £ A s H'_®A r H + . (2.11) 

r,s=0 

This space of multi-particle states is called the Fock space. 

There is a subset of states in this Fock space cosnisting of wedge products 
of single particle states. ( A general state is a linear combination of such 
wedge products.) Suppose, W is a subspace which doesnt differ 'too much' 
from H- . That is, the intersection of W with H+ and the intersection 
of W 1 - with 7i_ are both finite dimensional. Then we can form a state in 
T in which all states of W are occupied and all those in W ± are empty. 
Suppose e,i is a basis in W fl H + and and fj a basis in W 1 - fl 7i_ . 
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Then the state in T corresponding to occupying W will contain particles 
in WnH + and holes in W ± HH- : 

f 1 ' A/ 2 ' A--- Aei Ae 2 A e 3 A---. (2.12) 

(Here, /■?' is the dual basis in [W 1 UH-]' ). Moreover this is independent 
(upto multiplication by a complex number) of the choice of basis. Thus to 
each such subspace W C 7i there is a state (upto scalar multiple) in the 
Fock space JF . In other words we have an embedding of the Grassmannian 
into the Projective space VJ 7 ) . This is the infinite dimensionl version of 



the Pliicker embedding familiar from algebraic geometry. [f25| , |2f 

The Grassmannian describes this subset of states in Dirac's theory. The 
condition that the operator M = \& — e corresponding to a subspace W 
be Hilbert-Schmidt is precisely what is required for the above construction 
to go through. (It is not necessary to require that the projection operator 
7r + : W — > 7i + be finite dimensional as we did above: the precise condtion 
for the construction to work is that it be Hilbert-Schmidt.) What is special 
about this family of states is that they are coherent states: they minimize the 
uncertainty in the physical observables. Hence they are the states that have 
a sensible classical limit. The set of such coherent states is the classical phase 
space of the theory. That is why it is sensible to choose the Grassmannian 
as the classical phase space of our theory of hadrons. 

Dirac's theory was originally meant to describe electrons. We should 
think of it now as applied to quarks. Since the quarks are also fermions, this 
is reasonable. But there is one impotant twist: each quark comes in N 
colors: thus each state in Ti can be occupied by iV quarks, not just 
one. Moreover, only states that are invariant under the action of the group 
SU(N) are allowed: these are the states that describe the hadrons. The 
classical limit referred to above is the large N limit. 

The operator $ has eigenvalue —1 on states that are completely 
filled and +1 on states that are empty. The variable M measures the 
deviation from the vacuum state, in which all the negative energy states are 
completely filled. A baryon is essentially a color invariant bound state of N 
quarks. Thus the Dirac Sea of quarks can also be thought of baryons: 

each state is filled either by iV quarks or by a baryon. An infinitesimal 
disturbance from this can be thought of either as the creation of a meson ( 
a quark-anti-quark pair) or as the promotion of a quark from the negative 
energy sea to a positive energy state. This is an infinitesimal change if the 
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number of colors N is large. The operator M thus has to be determined 
by a map W : 7i_ — > 7i + . The block-diagonal elements of M are zero in 
this infinitesimal limit. 

In addition to such small deviations from the vacuum, the theory also 
allows for topological solitons which are not connected to the vacuum by any 
continuos path. These are very important, as they describe the baryons. To 
see how they arise we need to understand the topological properties of the 
Grassmannian. This is the subject of the next section. 

2.3 Renormalized Dimension of a Subspace 

Recall that $ = e+M has eigenvalues ±1 . Hence in some formal sense the 
trace of $ is an even integer: the difference between the number of positive 
eigenvalues and the number of negative eigenvalues. But of course this trace 
is not convergent. Even if we subtract the contribution of the vacuum, it 
is divergent: the trM can diverge, only the trace of its square needs to 
be convergent in general, since we only required it to be Hilbert-Schmidt. 
However the trace of M is conditionally convergent: we can define it to be 
the trace with respect to the above splitting into submatrices: 



This trace exists since a and d are trace-class matrices. It can be shown 
to be an even integer, which is invariant under the continuous deformations 
of the operator M . 

The interpretation of $ in terms of subspaces will help us understand 
the meaning of this integer. Imagine that we take a state from 7i + and add 
it to to get a subspace W : then W has dimension one more than 

7Y_ . This should change the trace of $ by —2 : one of the eigenvalues of 
$ has changed from +1 to —1. Since trM is essentially the difference 
between the traces of $ and e , we have tr e M = —2 . Thus 



is the 'renormalized dimension' of the subspace it describes: the difference 
between its dimension and the dimension of the standard subspace 7i_ . 

In fact this is the only topological invariant of M : any two oper- 
ators with the same conditional trace can be connected to each other by 



tr e M = tra + trd. 



(2.13) 



B = — trM 
2 



(2.14) 
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a continuous path. Thus the restricted Grassmannian is a union of con- 
nected components labelled by an integer, which is called the 'virtual rank' 
or 'renormalized dimension' of that component. Since it is invariant under 
all continuous deformations, in particular it will be invariant under time 
evolution. This is true for any reasonable definition of time evolution, inde- 
pendent of the choice of hamiltonian. Thus the renormalized dimension is a 
'topologically conserved quantity'. 

Each connected component of the Grassmannian by itself is an infinite 
dimensional manifold. Although the different components are the same (dif- 
feomorphic) as manifolds, the component of rank zero has a special role as it 
contains the 'vacuum' M — . An example with renormalized dimension 
one is a 'factorizable' operator, 

M = -2^®^ f . (2.15) 

The constraints on M are satisfied if ip is a positive energy state of length 
one: 

e^ = ^, |M| 2 = 1. (2.16) 

In the interpretation in terms of the Dirac theory, we have filled a positive 
energy state ip in addition to all the negative energy states. Hence this 
state contains an additional fermion occupying the state ip . 

More generally, let ip a be an orthonormal system of states satisfying 

eij = e a ij, <ij a ,ij b >=5l (2.17) 

Then 

M = -2j2^a®tpi (2.18) 

a 

is a solution of the constraints with renormalized dimension J2 a e a ■ 111 this 
configuration, we have filled a certain number if positive energy states ( ip a 
with positive fi a ) while creating holes in some others, ip a , with negative 
li a . The renormalized dimension is just the difference between the number 
of occupied positive energy states and the number of holes in the negative 
energy states. 

Thus we see that the renormalized dimension is just the fermion number 
if we interpret the Grassmannian in terms of the Dirac Sea. We will see 
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that the separable configuration is related to the valence quark model of the 
baryon. 

As noted earlier, a > and d < ; thus a contributes negatively 
to the baryon number and d positively. The only configurations that have 
6 = are the separable ones above. To see this note that, 6 = implies 
that a and d are proportional to projection operators. In order to 
be trace class, these have to be finite rank projections: corresponding to a 
certain number of baryons and anti-baryons. 

The valence parton approximation of the parton model corresponds to 
the separable ansatz in our soliton model. In that approximation, 4>a{p) is 
zero: there are no antiquarks in the proton within the valence approximation. 



2.4 Some Submanifolds of the Grassmannian 



In principle a baryon can be in any of the configurations of renormalized di- 
mension one: the ground state (proton) will be the one of least energy. Once 
we have determined the energy function (see below), the proton structure 
functions are determined by minimizing this energy over all possible config- 
urations. Later on, we will describe a method to do just that numerically: 
the steepest descent method. However this is a computationally intensive 
and slow method. A much faster method will be to minimize the energy 
over some submanifold of the Grassmannian, which is chosen so that the 
minimum in this submanifold is close to the true minimum, yet is easier to 



Physical intuition plays an important role in the choice of the ansatz 
of configurations explored this way. Such restricted phase spaces play an 
important role in our conceptual understanding as well: we will see that the 
theory restricted to rank one configurations of the type 



is just the valence parton model: we will be able to derive the parton model 
from QHD this way. 

Note that a change of phase ip — > e ld ip does not affect M . Indeed 
the set of rank one configurations is a submanifold of the Grassmannian 
diffeomorphic to V(H + ) . If we restrict the dynamics of our theory to this 
submanifold, we get an approximate theory which we will show is equivalent 
to the valence parton model. 



find. 





(2.19) 
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Allowing for a slightly larger set of configurations will give us the parton 
model with Sea quarks and anti-quarks. We will describe here the subman- 
ifold of the Grassmannian that descibes this approximation to our theory. 
We generalize the above rank one ansatz to a finite rank ansatz 

M = (2.20) 

Here, the ip a are a finite number r of eigenstates of e which are 
orthonormal: 

e^ a = e a ij a , ^ = 5% (2.21) 

The constraints 

M \r. (e + M) 2 = l, --tr e M = l, trM 2 < oo (2.22) 

2 

become the constraints on the r x r matrix ji : 

H = n\ [e + /i] 2 = l, ~tr^ = l. (2.23) 



Here, 



(ex 
e 2 ••■ 



is a diagonal r x r matrix, the restriction of 



V • • • e r J 

e to the finite dimensional subspace spanned by the ip a . The condition 
trM 2 < oo is of course automatic since M is now an operator of finite 
rank. 

Thus, given a set of vectors ip a satisfying the above conditions, we get a 
point in the Grassmannian of renormalized dimension one if fi itself belongs 
to a finite dimensional Grassmannian. If there are r + vectors ip a with 
positive momentum and r_ with negative momentum, we have 

tr[e + fi] = r + - r_ - 2. (2.24) 

The simplest solution is, r + = l,r_ = 0, so that fi is just a number: it 
then has to be —2 . This is the rank one solution desribed earlier, which 
leads to the valence parton model. 

If r = r + + r_ is two, there is no solution to the above requirements 
(more precisely all solutions reduce to the rank one solution). The next 
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simplest possibility is of rank three, with r + = 2, r_ = 1 . Then each 
solution to the constraints on // determines a one dimensional subspace of 
C 3 ; i.e., a point in CP 2 . In other words, for each such /x , there is a unit 



vector ( = 




E C 3 such that 



e ~ + /i = -l + 2C®C T , IICII =1- 



(2.25) 



Now, there is a C/(l) x C/(2) "gauge symmetry" in the problem: we can 
rotate the two positive energy vectors i/j and tp + into each other and 
change the phase of the negative energy vector ip_ without changing M , 
provided we make the corresponding changes in /i as well. This £7(1) x £7(2) 
action changes the phase of and rotates the other two components among 
each other. This freedom can be used to choose 



/ C- 



C->o. 



(2.26) 



In summary the submanifold of rank three configurations of renormalized 
dimension one is given by 

M = -2^o®'0o + 2 C-{C-['0-« ) '0- --0+<H>-0+] 

+ v /[i-C 2 ][^^l + ^+®^]} 



where ip-,ipo,i J + are three vectors in H satisfying 



(2.27) 



||^_|| 2 = ||^o|| 2 =||^ + || 2 = l, *+ = 0. (2.28) 

Also, 

< C- < 1- (2.29) 

The special case (- = reduces to the rank one ansatz. 

The physical meaning is clear if we consider the negative momentum 
components of M . The anti-quark distribution function is 

Up) = M aa (- P , - P ) = e\i>- a (p)\ 2 . (2.30) 
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Thus C 2 , is the probability of finding an anti-quark inside the baryon. 
Also, C-li'-ai— p)\ 2 is the distribution function of the anti-quark of flavor 
a . In the same way I^OaG )! 2 + ( 2 -\i J +a(p)\ 2 is the distribution function 
for quarks. In a loose sense \ipo a {p)\ 2 is the valence quark wavefunction 
and |^ +a (p)| 2 is the Sea quark wavefunction. But the splitting of the quark 
distribution into a valence and a Sea distribution is rather arbitrary and has 
no real physical meaning. 

By using the ansatz above with a larger and larger rank we can get better 
and better approximations to the Grassmannians. But the expressions get 
quite complicated beyond the rank three ansatz. We will find that this 
rank three ansatz holds the key to understanding the anti-quark and sea 
quark disributions in the quark model. A variational ansatz based on these 
configurations will provide a derivation of the parton model with sea quarks 
and anti-quarks from quantum hadrondynamics. 



2.5 Integral Kernels 

We can express the operator M in terms of its integral kernel in position 
space0: 

MiP(x) = J M(x,y)iP(y)dy. (2.31) 

Alternately we can think in terms of its action on the momentum space 
wavef unctions: 



$(p) = J ?p(x)e- ipx dx, VO) = / Hp)e tpx ^ (2.32) 



and 



[M^](p) = J M(p,q)i>(q)^. (2.33) 
The two points of view are of course related by Fourier transformation: 

M(x, y) = J M(p, g ) e *-<w^. (2.34) 



We will often suppress the indices a, b etc. labelling the basis in C N2f 
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The operator e is diagonal in momentum space: 

$>(p) = sgn (p)^(p). (2.35) 
It can also be described in position space as a distribution: 

eil>{x) = [ e(x - y)i>(y)dy — —V ( -^—^{y)dy (2.36) 
J it J x — y 

In fact this is a well-known object in complex function theory: the Hilbert 
transform (except for a factor of i ). It relates the real and imaginary parts 
of analytic functions. 

The constraints on M become then, 

M(p,q) = M*(q,p) 

~ f ~ civ 

[ sgn (p) + sgn (q)]M(p,q) + / M(p,r)M(r,q)— = 

f\MM\'*g < oo 

or equivalently, position space, 

M(x,y) = M*(y,x) 
J[e(x,y)M(y,z) + M(x, y)e(y, z) + M(x, y)M(y, z)]dy = 

J \M(x,y)\ 2 dxdy < oo. 

It is useful to see what the tangent space at the origin is like. Since M 
is infinitesimally small, the constraint can be linearized: 

[ sgn (p) + sgn (q)}M(p, q) = 0. (2.37) 

Thus M(p, q) is only non-zero if p and q are of opposite signs. In fact, 
the case where p > 0, q < determines the opposite one because of the 
hermiticity condition on M . Thus the tangent space is just the space of 
square integrable functions of one positive variable and one negative variable 
M(p, q) . 



2.6 The Infinite Dimensional Unitary group 

Given any self-adjoint operator $ , g^g^ is also self-adjoint. If the 
transformation g is unitary, gg^ = g^g = 1 , it preserves the condition 
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$ 2 = 1 . Indeed the action of the unitary group on such operators is tran- 
sitive: any operator satisfying $t = $ and $ 2 = 1 can be taken to any 
other by a unitary transformation. In terms of the variable M , the unitary 
transformation is 

M»gMgi + g[e,g*]. (2.38) 

In the infinite dimensional case, we have required that M be Hilbert- 
Schmidt (H-S). This means that there is a corresponding condition on the 
family of allowed unitary transformations: we define the restricted unitary 
group, 

U(H,e) = {g\gg^ = l,{e,g]el 2 }. (2.39) 

It is straightforward to verify that the H-S condition on the commutator is 
preserved under the multiplication and inverse operations, since X 2 is an 
ideal in the algebra of bounded operators. 

Indeed the restricted Grassmannian is a homogenous space of this re- 
stricted unitary group: 

Gr{l-L, e) = U{H, e)/U{H+) x U(H-). (2.40) 

The Lie algebra of the restricted Unitary group is 

U(H, e) = {u\u = -u\ [e, u] G J 2 }. (2.41) 

The infinitesimal action on the variable M is: 

M ^[u,e + M]. (2.42) 

2.7 Poisson Structure 

We will regard the Grassmannian as the phase space of our dynamical sys- 
tem. The matrix elements of the operator M are then a complete set of 
observables. We will seek a set of Poisson brackets among these variables. 
Fortunately there is a natural choice: there is a unique choice that is invariant 
under the action of the restricted unitary group. In fact the Grassmannian 
can be viewed as the co-adjoint orbit of the (central extension of the ) unitary 
group. We can regard the Poisson structure as induced by the Kirillov sym- 
plectic form on this orbit. However, the physical arguments use the Poisson 
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brackets rather than the symplectic form so in thsi paper we will not say 



much aboutthe sympletic form. (See Refs. 0, 14] for an elaboration of this 
more geometric point of view.) 

If the Poisson brackets are invariant under the action of the group, the 
infinitesimal action, 

M -> [u, e + M) (2.43) 

would be a canonical transformation for any u satisfying 

u = —u\ [e, u\ E I 2 . (2.44) 

The natural choice of a function that generates this canonical transformation 
is of the form 

/„ = k tr uM (2.45) 

for some constant k . There is a technical problem: the trace may not 
exist. Now remember that under the splitting TC = H- © H + , 

a p\ ( B I 2 \ ( a b\ (l l I 2 



U ={-pl 5) E {l 2 B)> M =W JHl 2 *) (2 ' 46) 
so that the conditional trace 

tr e wM = traa + tv/3b - trfttf + tv5d (2.47) 

exists. Moreover its value is some imaginary number, so that the constant 
k should also be purely imaginary in order that 

f u — k ti e uM (2.48) 

be a real-valued function. 

Thus we postulate the Poisson brackets 

{fu,M} = [u,e + M}. (2.49) 

These imply, of course that 

{fu, fv} = k tr e v[u, e + M] = -f[ u , v ] + k tr t v[u, e]. (2.50) 

The commutation relations of the restricted unitary group are satisfied only 
upto a constant term: it is the central extension of the unitary group that 
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acts on the Grassmannian, not the group itself. This has been studied at 
great length in the book by Pressley and Segal [27] so we wont go too far in 
that direction. 

We can write the Poisson brackets also in terms of the integral kernel of 
the operator M in position space: 

k{M^x,y),M c d (z,u)} = 5 c b 5(y-z){e a d (x,u) + M«(x,u)] 

-5 a d 5(x-u)[et(z,y) + M^z,y)], 

or momentum space: 

k{MZ(p,q),M c d (r,s)} = 5 6 c 2tt% - r)[8 a d sgn (p - s) + M d °(p, s)] 

-8P-k8{s - p)[5 c b sgn (r-q) + M c d (r, q)}. 

The principle of invariance under the unitary group cannot determine the 
constant k in the Poisson brackets. In a sense it doesnt matter what value 
we choose for it, as long as the values of all the observables are also multiplied 
by k . However it will be convenient to choose the value k = ~ as we will 
see in the next subsection. 

Soon we will have to calculate the Poisson bracket of functions f{M) 
which is not linear in M . The derivative of such a function can be thought 
of as an operator valued function f'{M) of M : 

df = txf{M)dM = J f a \M)dM^p,q)^ (2.51) 

Then we can use the above Poisson brackets to get 

k{f,M} = [f(M),e + M] (2.52) 

where the l.h.s. involves the commutator, as operators, of f'(M) and 
e + M . We leave the proof as an exercise to the reader. 



2.8 Momentum 

The above Poisson brackets are invariant under translations. There must be 
a canonical transformation that implements this symmetry. Formally, this is 
given by chosing u to be the derivative operator with respect to position; 
or, i times the multiplication by momentum. (Strictly speaking this is not 
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in the Lie algebra of the unitary group since it is not bounded, but let us 
ignore this technicality for the moment.) Thus, momentum is 

P = ik J pM(p,p)^. (2.53) 

Now, imagine calculating this quantity for the separable configuration of 
renormalized dimension one: 

M = -2i)^i)\ IH| 2 = 1, ety = i>. (2.54) 
With the choice k = | we would have, 

,dp 



poo 

p 



r p\^p)\ 2 4- (2 - 55) 



This has a simple physical interpretation: P is just the expectation value 
of momentum is a state with wavefunction tp(p) ■ 

Note that P is always a positive function: the quadratic constraint 
becomes, for p — q , 

2sgn (p)M(p,p) = - J \M{p,r)\ 2 — (2.56) 

so that momentum can be written as 

P = \j\p\ \M(P,r)\ 2 ^>0. (2.57) 

This makes physical sense if we regard this as the null component of the 
momentum vector. 



2.9 Kinetic Energy 

For a free particle of mass fi , the time-like component of momentum 
(energy), p is related to the spatial component pi by the condition 

p\-p\ = H 2 . (2.58) 

If we introduce the variables 

P = Po+Pi, P-=Po~Pi (2.59) 
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this becomes 

(2p - p)p = /i 2 , pp„ = fi 2 . (2.60) 

Thus the kinetic energy is related to the null component of momentum 
through the dispersion relation: 

Po = kp+~), P- = ~. (2.61) 
2 p p 

Notice that in this point of view the sign of energy p and of momentum 
p are the same: this will prove to be convenient when studying the structure 
of the ground state of relativistic fermion theories. (See the Appendix A for 
more details.) 

Thus we will postulate the kinetic energy of our dynamical system on the 
Grassmannian to be 

,2 



By the same argument as for momentum we can see that 

* = l/lW + Rl|tffor)|"!!>0 (2-63) 

If we add to this an appropriate potential energy U we will get the hamil- 
tonian of the system. This potential energy must transform like p_ , in 
order that (p,P- + U) transform like the null components of momentum. 
Under Lorentz tranformations, 

p — > Xp, p_ — > A _1 p_. (2.64) 

The Poisson bracket be invariant under Lorentz transformations, 

M(p, q) » M x (p, q) = XM(Xp, Xq). (2.65) 

In position space, x — > X^x so that 

M{x, y) M x (x, y) = X^'M^x, X~ 1 y). (2.66) 

Lorentz invariance will constrain the form of the potential energy as well. 
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2.10 Hamiltonian 



The total hamiltonian will be a sum of kinetic energy K and a potential 
energy U . While the kinetic energy is best understood in momentum space, 
potential energy is best written in position space. The simplest choice of U 
will be a quadratic function of M(x, y) . (Anything simpler will lead to 
linear equations of motion.) Thus we postulate 

U[M] = MZ(x,y)M b a (y,x)v(x,y)dxdy (2.67) 

+ -^J M^(x,x)M^y,y)v 1 (x,y)dxdy. 

This will have the requisite invariance under the internal symmetry U (F) . 

Lorentz invariance requires that U transform like p- or, equivalently, 
like - or x . Thus, 

v(\~ l x, \~ l y) = X~ 1 v(x, y). (2.68) 

Moreover, they can only depend on the difference x — y due to translation 
invariance. Thus 

1 1 

vi(x,y) = -a^x -y\, v 2 (x, y) = -a 2 \x - y\ (2.69) 

for some pair of constants a±, a 2 ■ 

Thus the hamiltonian of our theory is 

E[M\ = -\J[P+JW(P,P)^ 

+y J MZ(x,y)M b a (y,x)\x-y\dxdy 

+^ J MS(x,x)M b b (y,y)\x-y\dxdy. 

The constants m, ai,a 2 have to fixed later based on experimental data. 
They are the coupling constants of our theory. 

We have already shown that the kinetic energy is positive. The potential 
energy is manifestly positive when a±,a 2 > the way we have written 
it in position space. Thus the minimum value of E is zero and it is 
attained at the point M — .( This is in fact why we choose to parametrize 
our system by the variable M and not $ .) This of course lies in 
the connected component with renormalized dimension zero. In the other 
components,energy will have a minimum again, but it is not zero. In fact we 
will spend much time estimating this ground state energy. 



34 



2.11 Singular Integrals 



It will be convenient for later purposes to express the hamiltonian in mo- 
mentum space variables. This will require the use of some singular integrals. 



The basic singular integral is the Cauchy Principal value [28 



S(p) dp 



V / -i-^— = lim 

J a p — q 2n e^0+ 



q ~ e , f b 1 Sip) dp 
q+e \ p - q 2n 



+ 



(2.70) 



This exists whenever / is Holder continuous [] of order greater than zero at 
the point p = q . (Of course we are assuming that a < q < b ) . The basic 
idea of the Principal value integral is to cut-off the integral by removing a 
small interval of width 2e located symmetrically about the singular point 
and then take the limit as e — ► . (If the function vanishes at the point 
p = q the limit might exist even if not taken symmetrically: then we have 
an 'improper integral' rather than a 'singular integral'.) 

If we had cut-off the integral asymmetrically, the contribution of the 
region close to the singularity would have been 



Siq) 



+ 



1 dp 1 , . €\ 



—/(g) log A (2.71) 



p — q2ix 2i\ €\ 



This can take any value as we let t\ and t\ to go to zero. Requiring 
that ei = €2 removes this ambiguity. Indeed there are many other rules 
that could have been chosen: the particular one we choose must be justified 
by physical considerations: rather like the choice of boundary conditions in 
the solution of differential equations. In problems of interest to us there is 
a symmetry p — > —p (ultimately due to charge conjugation invariance) 
which selects out the Principal Value (or the Finite Part we will define soon) 
as the correct prescription. 

We already saw a use of the Cauchy principal value in the definition of 
the operator e : 

e^ x ) = l -V [ ^-dy. (2.72) 
7r J x — y 

The symmetric choice of regularization implicit in the Principal value is re- 
quired by the condition that e is hermitean. 

3 A function / is Holder continuous of order v (or / 6 C" ) if \im p ^ q ^^JZ^i q ^ 
exists for all q . 



35 



We will often have to deal with integrals that have a worse singularity: 

rb f(p) dp 



J a (p — 



(p-q) 2 2n 



(2.73) 



The 'Hadamard Finite Part' ('part finie') of such an integral is defined in 
terms of the Cauchy Principal value: 



TV 



f(p) dp 
(p — q) 2 2n 



V 



f(p) ~ /(g) dp 
(p — q) 2 2i\ 



(2.74) 



As long as / G C l , this will make sense: the vanishing of the numerator 
on the r.h.s. makes the integral exist as a Cauchy principal value. 

Although we wont need this yet, we note for completenes the definition 
of the finite part integral for a finite range of integration: 



TV 



b fjp) dp 
a (p — q) 2 2% 



V 



b f(p)-f(q) dp f(q) 
{p-q) 2 2tt 2tt 



a — q 



(2.75) 



We need the notion of a finite part integral because 



TV I \e ipx ^- = —\x-y\. 



2tt 



(2.76) 



The quantity on the right hand side is a Green's function of the Laplace 
operator on the real line: 



ill 

dx 2 2 



x-y\=S(x- y). 



(2.77) 



The choice of boundary conditions in this Green's function corresponds to 
the choice of the definition of the singular integral above. 
As an aside we note that 



TV 



g{p) dp 
p 2 2tt 



lim 



+ 



g(p) - #(0) dp 



2tt 



> 



(2.78) 



if 



g(p) > g(o) 

for all p . This is useful in checking that the hamiltonian is positive. 



(2.79) 
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2.12 Hamiltonian in Momentum Space 

Now, 

±jM;(x,v)Mt(y,x)\x-v\dxdy= (2.80) 

Jr(x-y)+ipx-iqy+ip'y-iq'x ^ X( ^y 

1 »>a/ \>>6/ .dpdqdr 



-TV J -MZ(p,q)M b a (q + r,p + r) 



Similarly, 



±J MZ(x,x)Mi(y, y )\x-y\dxdy (2.81) 
-TV / ^fe^.rt*!^* 

y e ir(x-y)+ipx-iqx+ip'y-iq'y^ x( ^y 

= -FP J ^K(p,p + r)M b b (p>,p> -r)^f 
Thus the hamiltonian becomes 

A J p 2n 



p 

~rP J ^M b %p, q )M b ( q + r,p + r)^ 



(2 



7T 



1 „~>a/ , \c*bt i i sdpdp'dr 



-^TV J -M:(p,p + r)M b b (p',p' 



(2k) 



In spite of the minus signs on the r.h.s., this is in fact a positive function 
of M . We already saw that the kinetic energy is positive whenever the 
constraint is satisfied. 

We can rewrite the second term in a more symmetric form, 



E[M\ = -± [\p+tL]M(p iP )± 
A J p 2ir 

n — -\ M b (n -I- — n —I— , 

2 1 (2tt) 3 



a\ f 1 ~„. r r,~ h , r r.dpdqdr 

~i TV l 7 M[ Xv ~ 2' q - 2 )MM + r p + 2 ] 
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-^p/-LM>,p + r) M t V,p'-r) Wr 

Now, 



(2vr; 



3 



satisfies 

/(r) = f*(-r) (2.82) 
so that the second term in the hamiltonian can be written as 

~tW> '<'>£■ (2 ' 83) 

Moreover, by a simple use of the Schwarz inequality, we see that 

Re/(r) < /(0). (2.84) 

Thus this term in the hamiltonian is positive. In the same way the last term 
can also be proved to be positive. 

2.13 The Equations of Motion 

From our Poisson brackets and the above hamiltonian we can derive the 
equations of motion. We can regard the derivative E'(M) — of the 
energy with respect to M as an operator valued function of M : 

dE = trE'(M)dM = J E' a b (M;q,p)dK(p,q)^. (2.85) 

Explicitly, 



dE = 
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Thus 



E' a b (M;q,p) = -hp + V-}2n5(p - q)8 b a 
4 p 

-YFVj-M^q + r^ + r^ 



(27T) 

-^l^TV - / M c (x/ v' + v-q) — 

2 °<^' {q-pyj m »W'P +P q >27r' 

Now the equations of motion implied by our Poisson brackets are, in 
operator notation, 

k^ = k{E(M),M} = [E'(M),e + M\. (2.86) 

(Here, k is the constant we fixed earlier to be | .) Substituting the 
above formula for E'(M) will give a rather complicated system of integral 
equations as our equations of motion. 

In particular, static solutions are given by 

[E'(M),e + M] = (2.87) 

The obvious solution to this equation is the vacuum, 

M = 0. (2.88) 

This is just the point where the hamiltonian is a minimum. Actually each 
connected component will have a minimum for the hamiltonian; the absolute 
minimum is the one with renormalized dimension zero. We will return to the 
study of the minimum in the components of non-zero renormalized dimension. 



2.14 Linear Approximation 

If M is infinitesimally close to the vacuum value M = we can linearize 
the equations of motion: 

-y-FP / h[M^p + s ,r + s)2ix5{r - q) sgn (q) 
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dvds 

-2irS(p - r) sgn (p)M fe a (r + s,q + s)}j—^ 



(p-r) 



2 2nS(r - q) sgn (q)M°(p',p' + r - p) 



;2n5(p — r) sgn (p)M^(p', p' + g — r) 



2tt (2tt) 



Simplifying, 

t dMZ(p,q) 

2 oft 
1 



—TV / -JAf? (p + s, q + s) sgn (g) - sgn (p)M t "(p + s, 9 + s)]— 



1 



2tt 



(p-9) 2 



(P - 

sgn (p)M c c (p', p' + q-p) 



sgn (g)M c c b',p' + g-p) 



dp' 
2^ 



Here, 



^(p) = i|p+^-] 

z p 



(2.89) 



Recall that the constraint on M becomes, in this linear approximation, 

[sgn(p)+ sgn(g)]M(p,g) = 0; (2.90) 

i.e., that p and q have opposite signs. We can assume that p > 0, q < 
since the opposite case is determined by the hermiticity condition 



Thus, 



1 dMg(p,q) 

2 dt 



M(p,q) = M*(q,p). 



--[K(p)-K(q)]M§(p,q) 

/I _ dg 
-Mg(p + s,q + s) — 



(2.91) 
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Now, under translation, M(x, y) — > M(x + a,y + a) and M(p, q) \— > 
e l(j '~ q ^ a M(p : q) . Thus the total momentum, which is a conserved quantity, 
is P = p — q . We can use as independent variables 

x B =p/P, P = p- q . (2.92) 
Clearly < x B < 1 and P > . For a stationary solution, 

M(p,q;t) = e^x(x B )- (2.93) 

We get thus 

-ux{x B ) = 

y + + pir^ 1 ^-* - ^* (2 - 94) 

Recalling that the mass fi' of the meson is given in light cone coompo- 
nents by, 

2cuP - P 2 = fx' 2 (2.95) 

we get 

V 2 x(x B ) = 

£ + jr^-Mxe) - ^FP f -^-5% - * /' x (»)* 
(1— xb) 7r Jo (x B — y) n Jo 

If we set 0:2 = this is exactly the equation that 't Hooft obtained for 
the meson spectrum of two dimensional QCD. 

't Hooft obtained this result by summing over an infinite class of Feynman 
diagrams. He showed that in the usual perturbaton series of two dimensional 
QCD, the diagrams of planar topology dominate in the large N limit. He 
then found a set of integral equations that describe the sum over these planar 
diagrams. 
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't Hooft's approach is limited by its origins in perturbations theory. It 
does not in fact reproduce the full large N limit of two dimensional QCD; 
only the linear approximation to the large N limit is obtained this way. 
The complete large N limit that we have constructed, by very different 
methods, is a nonlinear theory, reproduces the earlier theory of 't Hooft as 
just the linear approximation. 

Our theory of course contains much more information: it describes all 
the interactions among the mesons, since the hamiltonian is a non-linear 
function of the co-ordinates on the Grassmannian. Indeed it is not even a 
polynomial which means that there are interaction vertices at every order in 
perturbation theory. Some of these have been determined, but our analysis 
gives all the infinite number of vertices in terms of a single constant: the 
geometry of the Grassmannian fixes all these uniquely in terms of ati . The 
first such vertex was obtained by a summation of planar diagrams as well; but 
this approach will quickly get bogged dpwn by combinatorial complications. 
Before our work there was no indication that it was even possible to get a 
closed formula for all the interactions among the mesons in the large N 
limit of two dimensional QCD. 

The solution of this linear integral equation will give the structure func- 
tions of mesons. Experimental data exists on the structure functiosn of 
mesons as well. It is not from Deep Inelastic Scattering, but from Drell- 
Yan scattering. Preliminary results on the comparison of these structure 
functions with data are encouraging. However the focus of these lectures is 
the much more conceptually deep problem of obtaining the structure function 
of the proton. 

In fact our quantum hadrondynamics is not only a theory of mesons but 
also of baryons. We will show in the next chapter that there are topological 
solitons in our theory that describe baryons. 



2.15 Quantization of HadronDynamics 

We have so far described a classical theory, whose phase space is the Grass- 
mannian. The Poisson bracket relations for the M(x, y) play a role anal- 
ogous to the canonical commutation relations. There are two distinct ap- 
proaches to quatization of this theory: 

(i) Find a representation of the Lie algebra in terms of operators ona com- 



plex Hilbert space ( 'canonical quantization', pursued in Ref. [12|; and, (ii 
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Realize the wavefunctions as holomorphic sections of a line bundle on the 



Grassmannians ('geometrical quantization' pursued in Ref. [|TJ . 

In the first approach, we can use the complete classification of highest 
weight unitary representation of the (central extension of the ) infinite- 



dimensional unitary Lie algebra, due to Kac and Peterson [p9fl . These turn 
out to be described by Young tableaux just as in the finite dimensional case, 
except they may have infinite depth. The question arises which irreducible 
representation to choose. The quadratic constraint selects out one class of 
representations, whose Young Tableaux are rectangular, of width N and 
infinite depth. Such a representation can be realized as 

M(x,y) = ^:q i \x)q i (y): (2.96) 

where q, q^ satisfy the canonical anti-commutation relations. The represen- 
tation is carried by the vector space of color invariant states in the fermionic 
Fock space. Thus we recover exactly the light-cone gauged fixed version of 
two dimensional QCD with N colors as the quantum theory of hadron- 
dynamics. We have shown that two dimensional QCD and two dimensional 
hadrondynamics are exactly equivalent at all energies and for all values of 
the number of colors. We have already given the details of this argument in 
Ref.lO. 



In the geometric approach, we construct a line bundle of Chern number 
N on the Grassmannian. There is a connection on this line budle whose 
curvature is the natural symplectic form of the Grassmannian. The infinite 
dimensional Grassmannian is a complex manifold, just like the finite dimen- 
sional ones. Holomorphic sections of the line bundle exist when the Chern 
number is a positive integer. An inner product can be established on the 
vector space of these holomorphic sections by using ideas of Segal |!0| or the 
measure of Pickrell |31| . The observables M(x,y) can be represented as 
operators on this Hilbert space. Again the formulae we get are just the same 
as that of light-cone gauge fixed QCD. See |L4| for details. 

In addition it is sometimes convenient to quantize a restricted version 
of the theory where the phase space is some symplectic submanifold of the 
Grassmannian. This will give us some insight into the connection between 
our theory and the parton model. We will return to this in the next chapter. 

To summarize,we have a theory that doesnt just describe mesons to all 
orders in perturbation theory; we even have a theory of baryons. The topo- 
logically non-trivial solutions of the theory- the solitons- describe baryons. 
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In the next chapter we will study these solitons in more detail and obtain a 
theory of the structure functions of the baryon. 
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Chapter 3 



Baryons as Solitons and Their 
Structure Functions 

We saw that our phase space is a disconnected manifold, each connected 
component being labelled by an integer 

B = ~- tr e M. (3.1) 

Thus each connected component will have a minimum for the energy. These 
will describe stable static solutions of the equations of motion. They are 
called topological solitons. 

The minimum of the energy in each connected component with a given 
value of B desribes a stable particle, as the only particles with lower energy 
would be in a sector with a different value of B . It cannot decay into those 
states since B is a conserved quantity. 

We saw that hadron dynamics describes, in the linearized approximation, 
the spectrum of mesons of two dimensional QCD. What does a topological 
soliton describe? It is an old idea of Skyrme that the topological solitons of 
a theory of hadrons are baryons. This idea was revived in the mid eighties 
@, H and shown to be consistent with QCD in four dimensions. 

With this in mind, we should expect that the topological solitons of our 
two-dimensional theory describe baryons in the Deep Inelastic region. In 
the center of mass frame, a baryon will (due to Lorentz contraction) have a 
thin flat shape: one that can be described within a two dimensional theory. 
This offers us the possibility of solving one of the long-standing problems of 
particle physics: to explain the structure functions of hadrons as measured 
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in Deep Inelastic Scattering. By solving the static equations of motion, we 
will be able to determine the 'shape' of the soliton, which will then give us 
the dependence of the hadronic structure function on the Bjorken xb vari- 
able. The equations we need to solve are certain singular nonlinear integral 
equations. They are exactly solvable, so we will have to resort to numerical 
techniques. Indeed some ingenuity is needed even in the numerical part of 
this project: 'off the shelf methods do not work due to the singularities in 
the integrals. These singularities do not pose any new conceptual problems 
however: we are dealing with a finite quantum field theory. The basic def- 
initions of these singular integrals ( as dealt with an earlir section) go back 
to the nineteenth century work of Cauchy and Hadamard. 
The static equations of motion are 

[E'(M),e + M] = (3.2) 

where E'(M) is the derivative of the energy, computed above. It will be 
most convenient to view them as integral equations for M(p.q) , the integral 
kernel in momentum space. Of course we must solve these equations subject 
to the constraints 

[e,M] + + M 2 = 0; tr e M = -2. (3.3) 
The last condition picks out the sector with baryon number one. 



3.1 Quark Distribution Functions from QHD 

If we split the operator M into submatrices according to the splitting 

n = © n+ , 

M=(; t »), (3.4) 

the operator a : 7i_ — > 7i_ is positive while d : 7i + — > 7i + is negative. 
The baryon number B — — | tr e M is the sum of two terms: one from 
positive momenta and one from negative momenta: 

B = -\Hl [a aa (-p,-p) + daa(p,p)}^-- (3.5) 
2 „ Jo 2tt 
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The first term is always negative and the second always positive. Thus, 
<t>a(p) — \o-aa{— P, —p) can be thought of as the distribution function of anti- 
quarks in a hadron while 4>a{p) = ~jd aa {p,p) is the distribution function 
for quarks. The quantity 

<f>a(p) = Up) - Hp) (3.6) 

corresponds to what is usually called the valence parton distribution. Its 
integral over p and sum over a is equal to the baryon number. 

So the quark distribution function can be thought of as the sum of the 
'valence' and 'Sea' contributions: 



Mp)=<p V a(p)+<t> S a(p)- (3-7) 

The 'Sea' quark distribution function is thus just the same as the anti-quark 
distribution function: 



<p S a(p)=Up)- (3-8) 

We find this terminology a bit confusing and mention it only for purposes 
of comparison. The point is that there is no physical meaning to the splitting 
of the quark distribution function into a valence and a sea quark contribu- 
tion: only the sum has a measurable, physical significance. The anti-quark 
distribution does in fact have a physical significance. The above definition of 
a Sea quark distribution is thsu completely arbitrary and is not imposed on 
us by any symmetry such as charge conjugation invariance. In fact we will 
find later that there is another splitting which is more natural, where the 
Sea quarks have a wavefunctions are orthogonal to the valence quark wave- 
function. This at least respects the Pauli principle for quarks. Even that 
spliting is merely a matter of convenience of interpretation: all comparisons 
with experiment should be in terms of the measurable quantities <f> a (p) and 
<f>a{p) . 

Actually the distribution functions are usually thought of as functions of 
a dimensionless variable xb ■ If P is the total momentun (in the null 
direction) of the hadron, the momentum of the quark can be measured as a 
fraction of P . So we should actually write 

<Pa{x B ) =~Y^ d aa (x B P), S (x B ) = \ Yj d a a{-X B P)- (3.9) 
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These identifications of the parton distribution functions in terms of the 
diagional matrix elements of M a b(p,p') can also be seen in terms of the 
formula we derived in terms of QCD. 

The structure functions are supposed to vanish for Xb > 1 ■ Why should 
the r.h.s. of the above equation vanish for such x B ? For large N , the 

3 

structure functions will be of order e~ cN ^ at x B = 1 • So even though the 
wavefunction doesnt vanish at xb — 1 , it is exponentially small. A more 
proper analysis of the large N limit, taking into account semi-classical 
corrections will in fact reproduce wavefunctions that vanish at xb — 1 ■ 
The difference between these two versions of the large N limit is like the 
difference between the microcanonical ensemble and the canonical ensemble 
in statistical mechanics. 

3.2 Rank One Ansatz 

Before we set out to solve the above integral equations, it is useful to con- 
sider a simpler, variational approximation to them. A solution to the above 
constraints on M is the rank one (or separable) ansatz, 

M = -2V><g>'0 t - (3-10) 

We showed earlier that this is a solution to the constraints with baryon 
number one if ip is a positive energy wavefunction of length one: 

e^ = ^ EJ\Mp)\ 2 ^ = 1- (3.11) 

We are not claiming that the exact solution of the static equations of motion 
are of this form: only that this is a reasonable variational ansatz for it. The 
physical meaning of this ansatz will be made clear in the next section: it is 
equivalent to the valence quark model of the baryon in the large N limit. 

The technical advantage of this ansatz is that it solves the constraint; 
more precisely it reduces it to the much simpler constraint that ip is 
of length one. Geometrically, the set of states we are considering forms a 
projective space: changing ip by a complex number of modulus one does 
not change M . The separable ansatz is an embedding of the complex 
projective space CP(T-C + ) of the positive energy states into the component 
of the Grassmannian with renormalized dimension one. Instead of minimizing 
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the energy on the whole Grassmannian, in this approximation, we minimize 
it on this submanifold. 

The parton structure functions have a simple meaning within this approxi- 
mation: the structure function a (xe) is just the square of the wavefunction 
in momentum space: 

M*b) = \Mx B P)\ 2 . (3.12) 

The anti-quark structure function is just zero in this separable approxima- 
tion. 

The energy on this submanifold can be easily calculated: 

W) = E/f^ + fll^p)! 2 ^ (3-i3) 



Oil + a 2 



£ J \Mx)\ 2 \My)\ 2 \x-y\dxd y (3.14) 

ab 



where as usual, 



1>a(x) = I M P )e ipx ^. (3.15) 



Q Such a 'positive momentum' function has an analytic continuation into the 
upper half plane in the x variable. A simple choice would be 

Mx) = JxTiW (3 - 17) 

with the location of the double pole serving as a variational parameter. C a 
is constrained by the normalization condition that the length of ip is one. 
(A simple pole would have infinite kinetic energy, so we consider the next 
simplest possibility.) It is easy to see that there is a term in energy that 
scales like b^ 1 (coming from p \ip{p)\ 2 ) and the remaining terms scale 
like b . Thus there is a minimum. This encourages us to proceed to a more 
accurate determination of the wavefunction ip that minimizes the energy. 



: We will find it convenient to denote 

-2 



"i + ot 2 = g (3.16) 
to agree with previous papers. 
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The integral equations we have to solve are 

\[p + -}Mp) + rvtf [ v( P - q)4> a (q)^r = x Mp) 

2 p J 2ix 



p 

v(p) = -- 2 rr a ( P +q)M<i)p-- 

p A JO 2lX 



subject of course to ^>(p) = for p < and / °° \^>(p)\ 2 ^ = 1 • The 
second of the above integral equations is just the momentum space version 
of Poisson's equation of electrostatics: 

V"(x) = -r\x)^ a {x). (3.18) 

From Gauss' law, 

V{x) ~ -\x\ (3.19) 
2 

as \x\ — > oo : it is just the electrostatic potential of a unit charge located 
near the origin. Thus 

V(q) ~ -j 2 (3.20) 

as \q\ — > . This is why the integrals are singular: the finite part prescription 
is just a way of imposing the boundary condition that V(x) ~ ^\x\ at 
infinity in position space. We also note that for a solution centered at the 
origin, 

$(x)=tj)*(-x) (3.21) 
and V(x) = V(— x) . The boundary condition at the origin we impose is 

V(0) = (3.22) 

which translates to 

TV j V(q)^ = 0. (3.23) 
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3.3 Approximate Analytic Solution for /x = 



It will be useful to have an approximate analytical solution; even if it only 
works for a physically uninteresting region it will help us to validate our 
numerical method. 

Consider the singular integral (setting W(q) = q 2 V(q) ): 



Fvjv(q)4, a (p + q)^ = <pj 



W(q)MP + q)-W(0)Mp) dq 



2ir 



= 2 V 



W(q)[ip a (p + q)+Mp- ?)] - 2W(0)ij a (p) dq 



\ v f iy^ \-Mp + g) + MP - g) - 2^ a (p)] dq 



Q 



2tt 



1 

+ 2 



V 



1 



2[W(p) - W(0)] dq 



2tt 



i>aip) 



2ir 



The first integral is not singular any more. The integral in the square brackets 
in the last term is in fact zero by our boundary condition: 



7 



2[W{p) - W(0)} dq 
~q 2 27 



TV 



jv{q) 



dq 
2^ 



V(0) = 0. 



(3.24) 



Thus 



TV JvmUp + ^ = IJ fe(P + q) + ^ " q) ~ 24 ' M] §3.25) 

So far we havent made any approximations on this integral: just rewritten 
the singular integral in a better way. The main contribution to this integral 
ought to come from the neighbourhood of the point q = . So we should be 
able to approximate the quatity in the square brackes by its leading Taylor 
series approximation 



TV IV(q)4> a (p + q)^ 



q V{q) 2^ 



(3.26) 



Then our nonlinear singular integral equation reduces to an ordinary differ- 
ential equation: 



1 a 2 ~ 

A p 



(3.27) 
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where b is determined by the self-consistency constraint 



(3.2? 



This equation can be thought of as the nonrelativistic Schrdinger equa- 
tion for a particle of mass b^ 1 in a linear plus Coloumb potential. The 
only additional complications are that b is determined by the above self- 
consistency relation and also, we have the boundary condition if)(p) = for 
p < .The only non-linearity in the problem is in the equation determining 
b. 

Such non-relativistic models with linear plus Coloumb potentials have 
been used to describe mesons made of heavy quarks [j34| . But the physical 



origin of this Schrodinger equation is completely different in our case: we get 
this equation in momentum space and not position space. The eigenvalue 
A does not have the physical meaning of energy for us: energy is to be 
determined by substituting the solution into the formula for the hamiltonian. 
Moreover, we are studying the fully relativistic bound state problem. And of 
course we are studying baryons not mesons. Still it is encouraging that it is 
possible to reduce a fully relativistic bound state problem to a mathematical 
problem that is no more complicated than the non-relativistic case. 

The special case p = is particularly simple. In this extreme relativistic 
case, the solution is an Airy function: 

^(p) =C-Ai( P ~ 2A ) (3.29) 

C is fixed by normalization: Il^ip)] 2 ^ — 1 • The eigenvalue is fixed 
by the continuity of the wavefunction: since it must vanish for p < we 
require it to vanish at p = as well by continuity. Then 

A = --£i[2&( ai + a 2 )]* (3.30) 

where £ = —2.33811 is the zero of the Airy function closest to the origin. 
Now we fix the constant b by putting this back into the nonlinear self- 
consistency condition. We get 



2b 



J* Ai (O^P 



(27r)s[/« Ai^O^y- 
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The Analytic Approximation to the Wave Function with m=0 
0.8 




P g 



Figure 3.1: The approximate analytic solution for ji — 



A 



u i£ Ai (Qde 

V(27r)[/« Ai2(0de] 



-g ~ 0.847589^. 



Moreover 



^(p) = CAi 
We plot the solution so obtained below. 



(3.31) 



3.4 Numerical Solution with the Rank One 
Ansatz 

We now describe how to solve our integral equation numerically. First we 
replace the nonlinear integral into a recursion relation. 

±\p+(£-2\ a $ a+1 (p) + a 1 V rV s (p-q)^ s+1 (q)^- = (3.32) 
2 p Jo 2tt 



and 



V.(p) = -^/ °°*(P + 9W?)^- 



(3.33) 



We start with an initial guess Vo(p) , solve the linear integral equation to get 
a solution ipi(p) ■ Among all the solutions of this linear integral equation 
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Numerical and Analytic Wave Functions 




12 3 4 

Figure 3.2: Comparison of Numerical and Approximate Analytic Solutions 
for /j, — . 

we pick the one without a node; this happens to be the one with the smallest 
eigenvalue A . Then we calculate V\ (p) as above and then again solve the 
integral equation to get tp2 and so on till our iteration converges. 

We pick the nodeless eigenfunction at each stage since we expect our final 
answer for the ground state eigenfunction to have this property. It is just a 
coincidence that this happens to have the smallest eiganvalue A : in any 
case A does not have the meaning of energy. 

Of course to carry out this algorithm, we need to convert the above inte- 
gral equations into matrix equations. That is done by the quadrature method 
(described in the appendix C ) for singular integrals. The whole procedure 
can be implemented in Mathamatica quite well. 

We plot the numerical solution so obtained against the approximate an- 
alytical solution we got earlier for the special case fi = . 

The curves are qualitatively the same: the difference has to do with the 
approximations we had to make in order to get an analytic solution. The 
physicaly interesting case in fact does not have /i — : indeed we will see 
that the value of fi 2 is in fact negative. For example when m = , we 
have, u 2 = — — . It is of much interest to see if there is a variant of the 
analytic approximatin method above that applies to this more realistic case. 
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3.5 Quantization of the Rank One Ansatz 



We have considered elesewhere [O, HM the quantization of our hadron dy- 



namics, to recover QCD. It is also instructive to study the quantization of 
a simpler version of our theory, corresponding to the rank one ansatz. The 
quantization of this truncated theory gives a sort of approximation to QHD 
(hence QCD) which is of some interest in itself. We will see in the next 
chapter that this is just the valence parton model with interactions between 
the partons as predicted by QCD. 
Given a vector satisfying 

e ^ = -0, |H| 2 = 1 (3.34) 

we have an element of the Grassmannian of renormalized dimension one: 

M = -2i)®i)K (3.35) 

But this element remains the same if we change ip by a complex number 
of modulus one: ip — *> e t0 ip . Thus we have an embedding of the projective 
space V(7i+) into the Grassmannian. This is a sympletic embedding: the 
symplectic form on the projective space induced by this embedding is the 
same as the natural sympletic form on the projective space. Thus we can 
regard the separable operators as forming a 'reduced phase space' describing 
part of the degrees of freedom of our theory. 

A way to understand this is to consider ip{p) as a complex- valued 
observable on the effective phase space. In order to reproduce the Poisson 
brackets of M , these must satisfy 

{Mp),mp')} = o = {r a (p),r b ( P ')}, (3.36) 

{Mp),r b (p')} = -i2tt5(p - q)6 b a . (3.37) 
They satisfy the constraints, 

^(p) = 0; for p < 0, f°° m P )\ 2 ^- = 1. (3.38) 

JO Z7T 

We should regard the function ip(p) on momentum space to be the fun- 
dametal variable from which other observables such as 

i>(x)= r^{p)e^ (3.39) 

JO 2,71 
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can be obtained. 

Both the Poisson brackets and the constraint are simpler than the full 
theory, which is why we consider this case first. 

The Hamiltonian of our reduced dynamical system is obtained by putting 
the ansatz into E(M) : 

EiW = Eri + ;]W)lf (3-40) 

a Jo 2 p 2tt 
+ ai ^ a2 J2 J \ip a (x)\ 2 \i/j b (y)\ 2 \x - y\dxdy. 

ab 

These Poisson brackets and hamiltonian form a perfectly well-defined dy- 
namical system on its own right. We will now consider how to quantize this 
theory, and obtain the rules for the semi-classical approximation. 

We can quantize the theory by looking for operators satisfying canonical 
commutation relations: 

[k(p),k(p)\ = o = $ a (p),f(p% [kip)^\p)\ = ^8ip- P ')^Ai) 

As usual classical Poisson brackets go over to quantum commutation rela- 
tions: 

{A, B} — > —ih[A, B\. (3.42) 

In our case we will denote the parameter that measures the quantum correc- 
tion, analogous to H , by . In a minute we will see that this number N 
must actually be an integer. Th elimit iV — > oo is the classical limit. The 
constraint on the observables can be implemented by restricting attention to 
those states satisfying 



L 



V> (3-43) 
o In 

Now it is obvious that a representation for our commutation relations is 
provided by bosonic creation annihilation operators: 

[5»(p), hip')] = = [t(p), a ]b (p% [a a (p), aV)] = 2tt5( P - q)S b a . (3.44) 
with 

$a{x) = -^a ffi (x), ^\x) = ^^ a (x). (3.45) 



56 



Then the constraint becomes just the condition that we restrict to states 
containing N particles: 

f 00 a^ip)~a p (x)^ = N. (3.46) 
Jo 2n 

This is why iV must be an integer!. 

Now we know that we are dealing with a system of N bosons interacting 
with each other under the hamiltonian 

V 



+ ai + a2 N I a !a (x)a^ b (y)a b (y)a a (x)\x -y\dxdy. 



2 

The classical (or large iV ) limit we have been discussing so far is just the 
mean field approximation to this many-body problem. The semi-classical 
approximation will give us the leading corrections in the case of finite N . 
This basic insight is due to Witten, in a by now classic paper |J. 

What are these bosons? We will see in the next chapter that these bosons 
are just the valence quarks of the parton model, sripped of their color!. 
Quarks are of course fermions. However the wavefunction of the system 
must be totally anti-symmetric in the color indices beacuse of the condition 
that the state be invariant under SU(N) . Thus in the remaining indices 
the wavefunction must be symmetric: if we ignore color the valence quarks 
behave like bosons. 

Note that the momentum of the particles created by oJ(p) is always 
positive. Thus the total momentum 

P= f°° p a ^(p)a a (p)^- (3-48) 

JO 2ll 

is a positive operator. Indeed on a state containing N particles, 

|ai,pi; a 2 ,p 2 ; • • -a N ,p N >= a^{ Pl )a^{p 2 ) ■ ■■a^{p N )\0 > (3.49) 
P is just the sum of individual momenta, each of which is positive: 

P\a 1 ,p 1 ;a 2 ,p 2 r ■ -a N ,p N >= [pi+p 2 H PN]\ai, P i; a 2 ,p 2 ; ■ ■ ■ a N ,p N >(3.50) 

A general state will be, in this basis described a wavefunction 

10 >= f iP 1 7^4>( a i,Pi'i- ■ - a N,PN)\ai,Pi;- ■ -aN,PN > -(3.51) 
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An eigenstate of P with eigenvalue P will satisfy 

[Pi H Pn] <P(ai,Pi; ■ ■ -a N , p N ) = P0(ai,pi; • • -a N ,p N ). (3.52) 

Since each of the momenta pi are positive, it follows that they must each 
be less than the total momentum P : 

<Pi < P, fori = 1,...JV. (3.53) 

We will see in the next chapter that the Pi are the momenta of the 
valence partons: the proton is a simultaneous eignstate of H and P . 
We have just seen a very important point: our theory at finite iV (but 
within the approximation of the valence parton model) predicts that the 
wavefunction must vanish unless each of the parton momenta Pi are less 
than P . But there is no N in this inequality; so it must hold even in the 
large N limit!. 

We will be finding an approximate eigenstate of H , by a variational 
principle: a sort of mean field theory, the large N limit. The naive choice 
is a wavefunction which is a product of single particle wavefunctions. But 
such a naive version of mean field theory will violate the exact inequality 
we just established on momentum eigenstates. We should find our varia- 
tional approximation to the eigenstate of the hamiltonian, within the space 
of momentum eigenfunctions. 

Thus we assume that the wavefunction is approximated by a wavefunction 
that is just a product except for the constraint that the momenta add up to 
P : 

^(ai,pi; a 2 ,P2] • • -a N ,p N ) = 2n5(^2pi - P)^(ai,pi)^(a 2 ,p 2 ) • • -^(0^,^^03.54) 

i 

Thus the fraction of the momentum carried by each particle is less than one: 

i>(p) = 0, unless < ^ < 1. (3.55) 

That is how we recover the fact that the quark distribution function must 
vanish when the Bjorken variable is greater than one. This is a sort of semi- 
classical correction in the approximation. We must solve the classical 
equations of motion for ip{p) subject to the boundary condition that it 
vanish outside of the interval < p < P . This requires a modification of 
the variational ansatz but it it is possible to do that. 

This variant of mean field theory is rather like the micro-canonical en- 
semble. The naive mean field theory where only the expectation value of 
momentum is required to be P is like the canonical ensemble. 
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3.6 Valence Quark Distributions 



Now we come to the point of comparison of the calculated valence parton 
distribuion function againt the data. It is not in fact necessary to make a 
direct comparison with the data on Deep Inelastic Scattering. A generation 
of phenemenologists have extracted the valance parton structure functions 
from the data. More precisely they have assumed a parametric form such as 

(p(x B ) = Ax%(l - x B ) U2 [l + a x x B + a 2 x 2 B ) (3.56) 

for the parton distribution functions and fit to all known data points. Since 
the data consists of measurement of the DIS (Deep Inelastic Scattering) 
cross-section at different values of Q 2 this fit uses the convolution of this 
distribution function with the structure functions of the quarks calculated in 
perturbation theory. There are about 6 structure functions (corresponding 
to u,u,d,d,s quarks and the and gluons ) so altogether there are about 30 
parameters in addition to the fundamental parameters of perturbative QCD. 
Altogether there are about a thousand data points, coming from measure- 
ment of the cross-section for e — p , e — n , v — p ox v — n scattering 
at various energies. Thus the mere extraction of these distribution functions 
has itself become a subfield of particle physics |32], [19|, |33|]. The major groups 



seem to be in general agreement with each other, although a detailed ana- 
lyzis of the errors in their parameters is not yet available. There are recent 
attempts to estimate the systematic and statistical errors, but they seem 
incomplete. 

The comparison should be made with the weighted average over the u 
and d quark distributions,weighted so as to get an isospin invariant com- 
bination. This is because we have not yet done a collective variable variable 
quantization of the isospin degrees of freedom of the proton, so the distribu- 
tion function we are computing is the isospin invariant one. 

Another complication in the comparison is that not all the momentum in 
the baryon is carried by the valence quarks: it is known that only about half 
of the momentum is in the valence quarks, the rest being in the anti-quarks 
(or 'sea' quarks) or gluons. This affects our momentum sum rule. Essentially 
what it does is that the N in the sumrule gets replaced by N c s = y 
where / is the fraction of the momentum carried by the valence quarks. 
Until we do a complete calculation allowing for the contributions of the sea 
quarks and the gluons (i.e., without assuming the factorizable ansatz for M 
and without ignoring the scalar fields <pA ) we must treat / as a parameter 
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Wave Functions Versus x B 




Figure 3.3: Comparison of parton wavefunctions y/(f>{x) . The MRST 
distribution agrees with our prediction best when the fraction / of the 
momentum carried by the valence quarks is about 0.6 . 

and choose the value that gives the best fit. The good news is that it is the 
only parameter: m? is fixed to be zercQ, since the up and down quarks are 
known to have masses that are smal compared to the QCD scale ati . ( 
m u ~ 5 MeV and m d ~ 10 Mev while a x ~ 100 MeV .) 

The only other parameter of the theory is g 2 ; but since the distri- 
bution functions are dimensionless function of a dimensionless variable, the 
dependence on a\ cancels out. 

In comparison with experimental data we should keep in mind that we are 
in fact ignoring some of the constituents of the baryon: the valence quarks 
are knwon to carry only about half of the momentum of the baryon. The rest 
is in the sea quarks and the gluons. Thus we introduce a parameter / , the 
fraction of the momentum carried by the valence quarks. We find that for a 
value of / = 0.6 the wavenfunction predicted by us agrees quite well with 
experiment. (The value of being about a half is consistent with other ways 
of looking at the situation.) Thus we have solved the problem of deriving 
the quark structure functions of the baryon from QCD. 

It still remains to calculate the anti-quark and gluon structure functions. 

2 This by the way means that fi 2 = m? - ^- is negative. There is no contradiction 
here, since the quark is not a true particle: the mass of the baryon is still predicted to be 
real. 
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We will turn to that in later publications. In a later chapter we will give 
some ideas that help in solving this problem. 

In the next chapter we will see how to reconcile the soliton model with 
the parton model. 
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Chapter 4 

The Parton Model 



This chapter is based on the paper |L5] which in turn was an expansion of 



some short statements made in Ref. 12 



We have been describing the idea that the baryon is a soliton, essentially 
made up of an infinite number of mesons. But this seems to contradict the 
simple and successful picture of the baryon as a bound state of quarks. At 
the time the Skyrme model was revived, it was shown that the static (low 
energy properties) of the baryon that are so well explained by the quark 
model (magnetic moments, mass differences, flavor multiplet structure) can 
all be rederived in the soliton model 0. But the only explanation for the 
structure of the proton in deep inelastic scattering still appeared to be the 
quark-parton model. In these lectures, we are describing how to explain the 
structure functions within the soliton model as well. In fact we will be able 
to go beyond the simple minded parton model. 

It is an old idea in the study of the structure of the proton that the 
constituents-the partons- can be thought of as free particles. This then raises 
the question of how they bind to form the proton in the first place. It is the 
binding mechanism that determines the wavefunction of the partons and 
hence the structure functions of the hadron. In this section we will present 
a model of interacting partons. Perhaps surprsingly, it will turn out to be 
equivalent to the soliton model in the factorizable ansatz (and in the large 
N limit). Thus the parton model is merely an approximation to the soliton 
model. Later, by a deeper analysis of the soliton model, we will derive the 
sea quark and gluon distributions as well. 

It is at the moment impossible to derive the particle spectrum of four 
dimensional QCD directly. There are many attempts to do compute the 
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spectrum by direct numerical simulation of QCD, but they have to yet to 
surmount many obstructions. For example, the asymptotic freedom of QCD 
implies that the energy of all physical states are exponentially small compared 
to the cut-off. Thus extra-ordinary accuracy is required of all direct numerical 
computations of the energy, a problem that is only exacerbated by the large 
( in principle infinite) number of degrees of freedom in the system. However, 
it may be possible to understand the spectrum of two dimensional QCD 
by diagonalizing its hamiltonian numerically. Being a finite (rather than 
asymptotically free) theory numerical calculations are much more reliable. 



4.1 Part on Wave functions 



We will follow largely the discussion in |15j The valence parton model as- 
sumes that the proton is made of iV partons (quarks) which are fermions 
transforming under the fundamental representation of color. The idea of 
color of course was not present in the original version of this model, due to 
Bjorken and Feynman. Also it was not known then that the partons were 
spin half particles and hence fermions. But these were straightened out soon 
with the identification of the valence partons with the quarks, footnotelt is 
known that only about half of the momentum of the baryon is accounted for 
by these valence partons: the rest must be carried by the 'sea quarks' and 
by the gluons. 

We will work think of the parton wavefunctions as functions of the null 
momentum p = Po + Pi ■ In addition, the partons carries the quantum 
numbers of flavor and spin (together denoted by a ) and color a . Thus 
the wavefunction of a single parton will be ip{a, a,p) . Since p > for the 
null component of momomentum, we require that ip vanish for negative p 

A baryon is made of iV such partons so its wavefunction is a completely 
antisymmetric function ip(ai, azi,pi, a 2 , QZ2,P2', '"5 a Ni ^NiPn) ■ However, 
the baryon must be invariant under color: transform under the trivial repre- 
sentation of color SU(N) . This means that the wavefunction is completely 
antisymmetric in color alone: 

4>(ai, ai,pi; a 2 , a 2 ,P2] • • • ,a N ,a N ,p N ) = e aiiQ2i ... ajv ^(cti,pi; a 2 ,p 2 ] • • • ajv,^jv$4.1) 

In other words, the wavefunction is completely symmetric in spin, flavor and 
momentum quantum numbers. If we factor out color from the wavefunction, 
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the partons behave as if they are bosons. 



4.2 Hamiltonian 



The kinetic energy of a single parton is 




(4.2) 



Hence the kinetic energy of the partons making up a baryon is 



i-a N J ° i=l z Pi 



; o>n,pn)\ 



2 dp l ■ --dpN 
(2tt) n 



(4.3) 



In the early versions of the quark-parton model, they were treated as if 
they are free particles. That cannot be quite correct, since they must after 
all bind to form the baryon. The simplest idea would be to let them interact 
through a pairwise potential aiv(x) : 



o a i H / ~ Xj)\i>(ai,xu ■ ■ ■ ; a N ,x N )\ 2 dxi ■ ■ -dx N . (4.4) 



Now we must allow for a self-energy term as well in the hamiltonian. The 
point is that the effective mass of the partons may not be the same as their 
bare mass. Just as the electrons in a metal have a different effective mass 
from the free electron, the effective mass of the parton inside a baryon may 
be different from that of a 'bare' parton: there can be a finite renormalization 
of the mass. We can allow for this by replacing the /z 2 in the kinetic energy 
by an effective mass /i 2 . 

Since the partons (quarks) cannot be separated out to infinity, one could 
wonder what the 'bare' parton mass m means: it has no meaning as the 
mass of any particle that can exist in an asymptotic state. However, it is 
possible to give the bare parton mass a meaning in terms of high energy 
processes that do not involve asymptotic states. For example the weak in- 
teraction can be used as a probe of the quark masses: the weak decay rates 
of the quarks are sensitive to the bare masses. This gives a way to make 
sense of the bare parameter m . But for our purposes, what matters is the 
effective mass ji . 
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Thus the hamiltonian of the valence quark model is 



„.TZ„ Jo r-r 2 Vi 



N 



X too 

+ o a i H / H v ( x i - Xj)\i>{ai,xi] ■ • ■a N ,x N )\ 2 dx 1 ■ ■ -dx N . 

1 ai-a N Jo i+j 



The ground state of this hamiltonian is the baryon of the valence quark 
model. 

We havent yet decided what potential v(x) to use. There is ample 
evidence that the quark-quark potential is linear [[34]. In any case as in QHD, 
Lorentz invariance will lead to this choice n our lightcone co-ordinates. So 
we choose v(x) — . 

With this choice we get exactly the same hamiltonian we had in the last 
chapter, from the quantization of the rank one ansatz with playing the 
role of h . Thus the rank one approximation to the topological soliton 
model is equivalent to this interacting valence parton model. 



4.3 Hartree Ansatz 

The ground state of a many boson system can often be described by mean 
field theory: each boson moves in the field created by all the others. More- 
over, all the bosons can be assumed to occupy the same single particle state 
in this ground state. After the color is factored out, the partons in our model 
behave just like bosons. Hence we should be able to simplify the problem by 
making this mean field approximation. More specifically, 

N 

V>(ai,pi; ■ • -a>N, Pn) = ^^C^Pi - P) W_^{a-i,Pi)- (4.5) 

i i=l 

Here P is the total momentum of the system. Since p > we must 
impose 

ip(p) = (4.6) 

for p < . It then follows that p is less than the total momentum P . 
The wavefunction satisfies the normalization condition 
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and the momentum sum rule: 



N / P\^{P)\ 2 ^- = P. (4.* 



, P\^(P)\ 2 Tr- 
io Ztc 

The energy is calculated by putting this ansatz into the earlier formula: 

-f I v(x - y)j2\^(^ x )\ 2 J2\^( b ^y)\ 2dxd y- 

1 J -°° a b 

At this point we notice that we have exactly the same problem as in the 
last chapter. The energy of the rank one ansatz for M(p, q) is exactly the 
same as that of the parton model in the Hartree approximation. Thus the 
mean field theory of the valence parton model is our classical theory on the 
projective space CP{H + ) .. This interpretation of the rank one ansatz in 



terms of the parton model was already noted in Ref. [[[4 . 

What would be the leading deviation away from the large iV or classical 
limit? We argued in the last chapter in exactly the same theory (but viewed 
as the quantization of the rank one ansatz) that it amounts to the ansatz 

$(Pi,---Pn) = 2n5( Pl +p2---pN - P)J[4>(pi). (4.9) 

i 

There are now some correlations between the quarks: each momentum is less 
than P and also can be eliminated in terms of the others. The latter 
is a small effect; but the inequalities on the momenta, < Pi < P cannot 
be ignored. 
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Chapter 5 

Beyond the Valence Parton 
Model 

We saw that the separable ansatz in the soliton theory corresponds to the 
valence parton approximation. Now we see how to go beyond that and get 
the true minimum of the energy. We can improve on the valence parton 
model by adding Sea quarks and anti-quarks; this requires us to invent an 
improved variational ansatz. In another direction we can minimize the energy 
in the Grassmannian by numerical methods. We will describe ideas in both 
directions. 

5.1 A Co-ordinate System 

Recall that the equations to be solved are 

[E\M),e + M} = 0, (5.1) 
along with the constraints: 

[e + Mf = 1, tr e M = -2. (5.2) 

Define 

61 = 6-2^®^, (5.3) 

where ip is the minimum in the subset of separable kernels; or at least 
a good approximation to it. Rather than parametrize the soliton by the 
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deviation from the vacuum e it makes more sense to use the deviation from 
our approximate solution e\ which is in the same connected component. 
So we define a new variable M\ by 



e + M = e 1 + M 1 . (5.4) 

Then the constraints on M\ are 

[ei, Mi]+ + Ml = 0, tr e Mi = 0. (5.5) 

The trace of Mi is zero since the baryon number is now already carried by 
ei : 

- \ tr[d - c] = 1. (5.6) 

Let W be the subspace where e\ takes eigenvalue —1 : 

W = H-® < i) > (5.7) 

where < ip > denotes the linear span of ip . We have an orthogonal 
splitting of Ji according to this splitting. We can now represent M\ in 
terms of a co-ordinate system centered at the point e\ : 

m - ofli + zzT 1 -! [i + zzr i z \ , , 

M t - -i y zt[1 + zz ^_ 1 zt[1 + zzt ]-i z ) ^) 

Here, Z : W L ^ W is an arbitrary Hilbert-Schmidt operator. It is 
straighforward to check that this is a solution to the constraint equations. 
Moreover, in some finite neighborhood of e\ , all points of the Grassmannian 
can be represented this way. Of course this co-ordinate system will break 
down if we go too far away from e\ : we are assuming that the true minimum 
lies close enough to the approximate solution to be in this co-ordinate chart. 

Given a fixed tp , the operator Z above provides a co-ordinate system 
on an neighborhood of the Grassmannian. This co-ordinate system only 
covers a part of the connected component with baryon number one. It takes 
a countably infinite number of such charts (corresponding to different choices 
of ip ) to cover the whole space of baryon number one configurations. 
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5.2 The Method of Steepest Descent 



We need to minimize the function E(M) subject to the constraints 

[e, M] + + M 2 = 0, -- trM = 1, trM 2 < oo. (5.9) 

A simple method for minimizing functions of several variables is steepest 
descent: we start at some initial point and move along the straight line 
opposite to the gradient of the function at that point a small distance. Then 
we recompute the gradient at the new point and repeat the procedure. If the 
topography of the constant energy surfaces is not too complicated we will 
eventualy arrive at the minimum. 

The main complication in our case is of course the constraint: the gradient 
vector E'(M) is not tangential to the Grassmannian so if we move opposite 
to it we will leave the constraint surface. Even if we project out the tangential 
component of the vector, we will still leave the surface if we move along the 
straightline in that direction. The proper geometric solution to this problem 
is to move a small distance along the geodesic at M tangential to the 
gradient vector. Then we will recompute the gradient and find the geodesic 
at the new point. This is feasible beacuse the Grassmannian is a homogenous 
manifold and we can find the geodesic on it easily using its high degree of 
symmetry. 

More explicitly, recall that the tangential projection of the gradient vector 

is 

T = j[e + M,[e + M,E'(M)]]. (5.10) 

The vector 

Y = ±[e + M,E'(M)] (5.11) 

is at right angles to the tangential part of the gradient: it is in fact obtained 
by multplying T by the complex structure of the Grassmannian. Just 
as in the case of the sphere, the geodesic is obtained by rotating the point 
about an axis orthogonal to the tangent vector. More explicitly, the geodesic 
starting at M tangential to T is 

7 (r) = e rY [e + M}e- TY - e. (5.12) 
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Thus the steepest descent algorithm is 

1. Choose an initial configuration Mq and small parameter r . 

2. Given the k th configuration M k , calculate the gradient E'(M k ) and 
Y k = \[e + M k ,E\M k )}. 

3. Set M k+ i = e rYk [e + M k ]e~ rYk — e and repeat the previous step. 

The value of r has to be chosen by some trial and error. Too small 
a value will produce changes in the configuration within the noise due to 
numerical errors. Too large a choice will not give a convergent sequence: 
we will bounce around all over the Grassmannian. But we found that in 
practice, a proper value of r can be found quickly. 

How do we choose the initial configuration Mq ? Since we believe 
that the valence approximation is good, we could use as the starting point a 
separable configuration, minimizing the energy within that subspace. 

5.3 Sea Quarks and Anti-quarks 

We saw that if we restrict the dynamics of our theory to the rank one ansatz, 
M = —2tp <S> ^ we get the valence parton model. To get a more general 
picture that includes Sea quarks and anti-quarks (but is still not the total 
picture) we must use the ansatz with rank three. In fact by going to ansatzes 
of larger and larger rank we can get better and better approximations to the 
baryon structure functions. 
The rank three ansatz is 

-V[i-Ci][^-®4 + ^+®^-]} 

where ip-,ipo,ip + are three vectors in 7i satisfying 

ei\)_ = e^po = ipo, = (5.13) 

|hM 2 = ||^oir = ||^ + || 2 = l, *+ = 0. (5.14) 

Moreover, 

< C 2 < I- (5.15) 
We saw that C- is the probability of finding an anti-quark inside the baryon. 
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As noted earlier, the anti-quark distribution function is just ( 2 \ip a (~ P) 1 2 
The baryon number is 



B = Ef {l ^°° )|2 + C-[l^(^)| 2 - ^-a(-p)| 2 ]}^. (5.16) 
- Jo 271 



P = T, rp{WM\ 2 \ + e[\^a(- P )\ 2 + \^a(p)\ 2 }}^ (5-17) 



These confirm the interpretation of ipo as the valence quark wavefunction 
and i/j+(p) as the Sea quark wavefunction. 
Also, the kinetic energy is 



K = T.[~ 2 \P+ ^-}{\^a{p)\ 2 \ + CU-a{-p)? + \^a(p)\ 2 }}^ (5.18) 



The potential energy is more complicated. Recall that the potential en- 
ergy is simpler in position space while the kinetic energy is simpler in mo- 
mentum space. Our wavefunctions now depend on a discrete variable in 
addition to momentum, and there is a unitary transformation in these dis- 
crete variables that is the counterpart to Fourier transformation. This is the 
transformation to a basis in which fj, is diagonal; the wavefunctions will no 
longer be eigenstates of e . 



The total momentum is, similarly, 




(5.19) 



where 



^2 



(5.20) 



(5.21) 



The potential energy is then, 



U 
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e\^a(x)\ 2 V U (x) + e-\i>2a{x)\ 2 V 22 {x)}dx 
+ Y 2Re E / {C-^Oa(xWla^)Vl (x)- 
C-1p0a(x)l^2a( X ) V 20(x) - C-^la{x)^l a {x)V 2 l{x) Jcfo. 

The mean fields are determined by solving the differential equations 

V^(x)=^ aa (x)^ a (x). (5.22) 
with the boundary conditions 

V aj3 {x) -> 5 Q/ 3-|z|, for |x| -> oo. (5.23) 

Equivalently 

Vaa(0) = C*(0) = (5.24) 

for the diagonal components, and, 

V a p(oo) = V^(oo) = for a + (3 (5.25) 

for the off diagonal components. 

By choosing appropriate variational ansatzes we can estimate the anti- 
quark content of the proton. Or we can derive integral equations for the func- 
tions ip±fl and solve them numerically. Both methods are being pursued. 
Initial results are encouraging. Detailed results will appear in a separate 
publication. 
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Appendix A 

A Null Co-ordinate System 



A.l Kinematics 

It will be convenient to use a co-ordinate system that combines the advan- 
tages of the null and Cartesian co-ordinate systems. 

If x° and x 1 are the usual Cartesian co-ordinates in Minkowski space, 
the metric is 

ds 2 = [dx ] 2 - [dx 1 ] 2 . (A.l) 

We define now 



x = x 



(A.2) 



so that 



ds 2 = dt[dt + 2dx]. (A. 3) 

The Minkowski metric is, 

V»»=( 1 1 J) (A.4) 

Thus the vector is time-like while ^ is null. The initial values of 

fields will be given on a surface of constant t , which is a null line. 

Momentum p = p^dx^ = podt + pdx is a 1 -form (co- vector); we 
will use the same letter to denote the momentum 1-form as well as its null 
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componenet, but it should be clear from the context which one we mean. To 
find its magnitude we must use the inverse of the above metric tensor: 

»T=(5 _\)> p 2 =p,p u r ] ^ = 2p p-p 2 . (A.5) 
Thus the mass shell condition becomes 

(A.6) 

Here, \i is the rest mass of the particle. 

We see now one of the main technical advantages of using the null-time 
co-ordinate system: energy and momentum have the same sign. In the usual 
space-time co-ordinates, p Q = ±-\/[Pi + A* 2 ] an d therefore no such simple 
relationship exists. In the Dirac theory of fermions the states of negative 
energy are occupied; this becomes merely the condtition that the negative 
momentum states be occupied. 



zr 

P+ — 

v 



A. 2 Dirac Matrices 

The Dirac matrices are best thought of as matrix-valued vectors, since they 
appear in the form 7 M V M in the action. Thus the Dirac algebra 7^7" + 
7^7^ — 2r) tiU becomes 

[ 7 f = 0, 7*7 + 77* = 2, 7 2 = -l- (A.7) 
We will choose the explicit representation 

*-(! «)• ■>-(! 7)- w 

In any representation of Dirac matrices there is a 'charge-conjugation' matrix 
C satisfying 

C-fC- 1 = ('f) T . (A.9) 

In the usual space-time formalism this matrix is often 7 itself, but that is 
a representation-dependent fact. In our representation, 

C=(\ J). (A.10) 
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Given the Dirac spinor q — (j^j > the conjugate spinor is 

q = q ] C={q\ q\ ) . (A.ll) 

A. 3 Free Fermions 

The Lagrangian of a free Dirac fermion becomes 

Ld = qY[-id^]q + mqq 

+2ql[-id t ]q 2 + q[{-id x )qi - q\{-id x )q 2 + m[q\qi + q\q 2 }. 

We see that the field qi has no dynamical degrees of freedom: it has no 
time derivative in the action. Hence it can be eliminated by its equation of 
motion: 

Tfl 

qi = -—ft (A.12) 
V 

where p = —id x . Putting this back into the action and changing variables 

X = 92 (A.13) 
gives us the effective action for the propagating field: 

1 TH? 

L x = x\-id t )x-x\[P+— ]X- (A.14) 

A. 4 Gauge Fields 

The lagrangian of two dimensional QCD is 

N 



L = N 



N r f r 

— / trF^F"" + J2 / q aa Hl ■ V + m a }q aa . 

We have added in the flavor indices a, 6 and the color indices i,j . 

The freedom of gauge transformations can be utilized partially to impose 
the null gauge condition A x — . Then, the fermionic part of the lagrangian 
becomes 

L D = ^Y[-i(d»6i + A^)]q aj + mq ai qai 
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+ql m (-td x )q lai - \x ]a \-id x )Xai + ^[x^Qiai + ?5"x«]- 
which becomes 

L x = X^(-id t ) X a ~ X ]a {^[P + y] - <^}x« (A.15) 

upon eliminating (ft . 

To this we must add the action of the Yang-Mills field itself, which looks 
quite simple in this co-ordinate system and gauge: 

Ly M = ^- tTF^rfrr = tr[d x A t } 2 . (A.16) 

The field A t does not propagate and can be eliminated. Thus the action 
of two-dimensional QCD can be written entirely in terms of the field x ■ 



A. 5 The Dirac Vacuum 

Let us return to the free fermion theory with Lagrangian 

j Tfl^ 

L = x\-idt)x-X^[f>+—]x (A.17) 

Upon quantization, the field x becomes an operator satisfying the fermionic 
anti-commutation relations 

lx(x), x(y)U = S(x - y), [ X (x), x (y)U = 0. (A.18) 

If there were only a finite number of operators, such canonical anti- 
commutation relations would have a unique representation. In the infinite 
dimensional case physical ideas have to be brought in to choose the right 
representation. Dirac showed that the correct choice is to assume that all 
the negative energy states are filled even in the vacuum. Since energy and 
momentum have the same sign in our co-ordinate system, this condition is 
easy to implement. We define the vacuum by 

x j {p)\0 >= forp < 0, x(p)|0 >= forp > 0. (A.19) 
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We have defined the Fourier transforms 



X(x) = J X{VY VX % (A.20) 

etc. Then the normal ordered product of a pair of operators is defined as 

■xKp)x{p')-=x\p)x{p') (A.21) 
unless both p and p' are negative, in which case it is 

: X t ip)x(p') ■ = -X{P')X\P)- (A.22) 

The point is that then the expectation value of normal ordered current op- 
erators are zero in the Dirac vacuum. Indeed, 

<0| X t (x)x(y)|0> = J '^e"*^ < 0|x f (p)x(<z)|0 > 

= f_^e^= l -[5{x-y) + e{x-y)]. 

Here, 

e(x - y) = ! sgn (p)e ip(x - y) ^- = —V—^— (A.23) 
J 2n m x — y 

is to be thought of as a distribution. It is (upto a factor of % ) the kernel of 
a well-known integral transform, the Hilbert transform. 

Thus we should regard the current operators of the fermionic theory as 
defined with the normal ordering. For example, the equation of motion of 
the A t will be, in the quantum theory, 

-%A> s (x) = ±a 1 :x* ai (x)xv(x):. (A.24) 

We can use this to eliminate the gauge field from the theory completely. Thus 
two dimensional QCD can be written as a theory of fermions interacting with 
each other through a Coulomb-like potential. The further analysis of this 
theory is carried out in the text, towards the end of the first chapter. 
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Appendix B 
Operator Ideals 



Here we give the basic definitions of the operator ideals we use in the text. 



A deeper discussion may be found in Ref. E7L 35, B6[; we give only a rough 



outline of the theory. It is the author's fervent hope that experts in functional 
analysis will not read this appendix. 



B.l Compact and Hilbert-Schimdt Operators 

The rank of an operator A : 7i — > Ti on a complex Hilbert space is the 
dimension of its range; i.e., the dimension of the subspace of all vectors that 
can be written as Au for some u G W . When the rank of A is finite, it 
can be thought of as a sort of 'rectangular matrix' with (possibly) an infinite 
number of columns but only a finite number of rows; at least there is a basis 
in which it has this form. 

For such an operator, we can define several measures of its size (norm). 
For example (the operator norm), 

\A\ = sup Mf. (B.l) 

u \\U\\ 

Another (the Hilber-Schmidt or H-S norm) is 

\A\ 2 = tx[AA^. (B.2) 

\A\\ is also the sum of the absolute squares of all the matrix elements in 
any basis. 
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The completion of the space of finite rank operators in the operator norm 
is the space of compact operators. In other words, a compact operator is one 
that can be approximated arbitrarily closely by finite rank operators, distance 
between operators being measured with the operator norm. If we instead 
complete in the H-S norm, we get the space of Hilbert-Schimdt operators. 

A bounded operator is one which has finite operator norm; i.e., 

II An 1 1 

|A| = sup^^ (B.3) 

uen \u\ 

exists. The set of bounded operators on Ti is an algebra B(H) . Not all 
bounded operators are compact; for example the identity is bounded yet not 
compact. 

A compact operator has the expansion, with fj, n > , 

A = X>n|^n >< n | (B.4) 

71=1 

the sum being either finite (when A is of finite rank ) or infinite (more 
generally). The numbers fi n are called singular values; if A is positive 
they are its eigenvalues. A and its adjoint A^ have the same singular 
values. 

Roughly speaking, the operator norm of A is the largest of its singular 
values, while the Hilbert-Schmidt norm is the sum of the squares of the 
singular values. The singular values of a compact operator form a sequence 
that converges to zero; in fact zero is the only limit point of the sequence. 
Thus, if the H-S norm is finite, the eigenvalues must be tend to zero. In fact 
all H-S operators are compact. 



B.2 Schatten Ideals 

Many other norms can be defined in terms of the singular values. For example 
the trace class operators are those for which the sum of singular values is 
convergent. This is stronger than the requirement that the diagonal elements 
in some basis form a summable sequence: being trace-class reqires a sort 
of absolute convergence. The product of two Hilbert-Schimdt operators in 
trace-class. 

More generally, for p > 1 , we have the class of operators for 

which the sum J] n v n converges. Of course X 2 is the space of H-S operators 



82 



and X\ that of trace class operators. The space of compact operators can 
be thought of as the limiting case X^ and that of finite rank rank operators 
as the opposite limit, X . We have the inclusions 

J C X 1 C X 2 ---X OQ C B. (B.5) 

It is very important for us that the X p are two-sided ideals in the algebra 
of bounded operators; i.e., that A G B and B G X p implies that both 
AB G X p and BA G X p . These are called the Schatten ideals. We are 
especially interested, of course, in the cases p — 1, 2 . 



B.3 The Restricted Unitary Group and its 
Grassmannian 

In the text we are interested in the case of a Hilbert space 7i with a 
given orthogonal splitting into two infinite dimensional orthogonal subspaces: 
H = H- © H + . Recall that the operator e is defined to have eigenvalues 
±1 on H± . The restricted Unitary group is the subset of all unitaries 
satisfying a convergence condition: 

EMW,e) = M = 9 ] 9 = 1; M e X 2 }. (B.6) 

It is vital for this definition to make sense that X 2 is an ideal; that is why 
the product of two elements still satisfies the convergence conditon. (Any 
unitary operator is bounded). If we split g G Ui(H,e) into submatrices 

according to the splitting g = ^ ^ where a : H- — > etc., the 

submatrices b, c are H-S. The matrices a, d may not be invertible in 
general, but they are Fredholm (see below) since g' 1 exists. 

The restricted Grassmannian is the set of all self-adjoint operators of 
square one whose distance ( in the H-S sense) from e is finite. 

Gn(H, e) = {$|$ f = $; $ 2 = 1, $ - e G J 2 }. (B.7) 

In the text we often use M = $ — e as the variable that describes a point 
in the Grassmannian. We showed that each such operator corresponds to a 
subspace of 7i ( the negative eigenspace of $ ) which is at a finite distance 
from the standard subspace 7i_ . 
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The fact that X 2 is an ideal ensures that the action of the restricted 
unitary group on this Grassmannian well-defined: 



$»g$g\ M^gMg^+g{e,g\ (B.8) 



B.4 Fredholm Index 



The material in this section is contained in Ref. |27j to which we refer for 
proofs and more precise statements. A bounded operator A is Fredholm if it 
is invertible modulo a compact operator; i.e., if there is a compact operator 
K such that A + K has an inverse. The set of Fredhom operators is 
a toplogical space, with the toplogy induced by the operator norm. It is 
disconnected, each connected component being labelled by an integer called 
the Fredholm index. 

To understand the Fredholm index, consider the kernel of an operator; 
i.e., the subspace of all u such that Au = . Even in the finite dimensional 
case the kernel can change change discontinuosly under a small change in A 
.; for example we may cross an eigenvalue. In the finite dimensional case, 
the dimensions of the kernels of A and A* are the same. In the case of 
Fredholm operators on infinite dimensional Hilbert spaces, the difference 

index (A) = dimker A — dimkerA^ (B.9) 

is always finite but need not vanish. It however does not change under 
continuos changes in A : it is constant in each connected component of 
the space of Freholm operators. In fact it is the only such function: any two 
Fredholm operators of the same index are connected by a continuos curve. 

Now g G Ui(7i,e) is of course Fredholm of index zero: it is invertible 
as is its adjoint. But the submatrices a, d ( where a : H- — ► 7i- etc.) 
defined above are not invertible in general. Yet they are invertible modulo 
some compact operators ( products such as bb* ), so they are Fredholm. 
Since g as a whole invertible, they must have opposite Fredholm indexes. 
The index of a (for example) is a topological invariant of g . The group 
U\ is the union of connected components labelled by this index, which can 
take any integer value. The connected component of the identity, of course, 
has index zero. 

We can understand the renormalized dimension of a subspace (called vir- 
tual rank in [PJ) in terms of the Fredhom index. Every self-adjoint operator 
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can be diagonalized; an operator $ € Gr\ can be brought to the standard 
form e by an element in U\ : $ = gegt . The index of a (g) is then 
the topological invariant associated to $ . 

Another point of view is in terms of the subspace W of H correspond- 
ing to $ . If $ is at a finite distance from e , W will not differ 'too 
much' from 7i_ . More precisely, the projection operator 7r + : W — > H.+ 
will be compact and 7r_ : — > 7i_ will be Fredholm. The index of 7r_ 
measures the 'difference in dimensions' between W and 7i_ ; this is the 
renormalized dimension of W . We saw in the text that this has the physical 
meaning of baryon number. 

It si crucial for all this that we allow only subspaces at a finite distance 
from the standard one in Gr\ . If we had allowed for all subspaces, the 
Grassmannian would have been contractible!. The convergence condition 
we must impose is required for the Poisson brackets to make sense: the 
symplectic form of the Grassmannian would not make sense otherwise. It is 
gratifying that as a consequence, we get a topological invariant which has 
the physical meaning of baryon number. 

Although the symplectic form makes sense on all of the phase space Gr\ , 
we should not expect the hamiltonian to make sense on all of it. There should 
be some dense domain in which the hamiltonian does make sense however. 
We leave such questions as challenges for the analyst who is interested in 
solving problems of relevance to physics. 
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Appendix C 

Quadrature of Singular 
Integrals 



In this appendix we summarize some ideas on the numerical methods that 
are used in the chapter on solitons. Some originality is needed even in this 
part of the problem. 



C.l Quadrature Formulas 

The equations we have are just too hard to be solved analytically. We have 
to resort to numerical methods. The basic idea is to convert the integral 
equation into a matrix equation by allowing p to take just a finite set of 
values: we must find a way to approximate the integral by a finite sum. Then 
we will solve the resulting nonlinear matrix equations by iteration. 

Quadrature is the approximation of integrals by finite sums. There are 
standard methods for quadrature, (method of moments) going back to the 
days of Gauss. But we have to modify these methods since our integrals 
are singular. The basic idea is well-known in the literature on quadrature at 



least for the case of the Cauchy principal value |37j. Our integrals are one 
step harder (Hadamard Finite Part) but the idea is the same. See S. Chan- 
drashekhar/s classic book [[$8[] for a clear discussion of numerical integration. 

Let p(x) be a continuous positive function on the close interval [a, b] 
. We are interested in evaluating integrals such as 

b 

f(x)p(x)dx (C.l) 
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by numerical approximations. Here f(x) is some continuous function. 

Let Xj, j = 1, • • • v be a set of points in the interval [a, b] . We expect 
a weighted sum such as 

»,/!•'•,) (C2) 

3 

to be good approximation for the integral, provided that (i) the number 
points v is large enough and (ii) the points X j £1X6 distributed roughly 
uniformly. 

Every function can be approximated by a polynomial within the interval; 
as the order v — 1 of the polynomials grows the approximation gets better. 
We can thus approximate the integral of f(x) by that of its polynomial 
approximation of order v — 1 . The weights Wj are determined (once the 
points Xj are chosen) by this requirement: the above formula should in 
particular be exact for polynomials of order v — 1 . 

This is the same as the condition 

rb 

/ x k p(x)dx = ^WjXj, for k = 0, 1, • --k - 1. (C.3) 

Ja ■ 

The left hand side are the moments of the distribution p(x)dx , and are 
assumed to be known. Then the above set of linear equations determine the 
weights in terms of Xj . 

If the function p(x)dx is not too rapidly varying, a simple choice such 
as equally spaced points Xj = a + (b — a)— should give reasonable 
approximation to the integral: certainly for polynomials upto order v — 1 
we will get the exact answer anyway. But it is possible to do better by 
choosing the v points Xj cleverly, as pointed out by Gauss: we can 
ensure that the answer is exact for polynomials of order 2v — 1 . But we 
wont be using this idea: we will just use equal spacing for the points, which 
turns out to be more convenient. This because our integrands involveterms 
such ip(p + r) , so it is convenient if the sum two points Xj + Xk is also a 
point at which we evaluate the integrand. The loss of precision in quadrature 
can be made up because the simplicity of equal spacing allows us to choose 
a larger number of points. 
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C.2 Singular Measures 



Now consider integral|37J TV J b f(x)^§ ■ The symbol TV indicates as 
before the Hadamard 'finite part' of the integral. For such singular integrals 
we can also find a numerical approximation as above. But the system of 
moments is no longer positive. This is because the integation measure is 
no longer positive: TV Jq f{x)^s can be negative even when f(x) is 
positive. The moments of the measure are given by 

TV \ x k ^- = f— - for k ± 1 (C.4) 
Jo x 2 k — 1 



and 



rb fjnr- 

TV x k ^r = \ogbioik = l. (C.5) 

JO X 2 

Note that the zeroth moment is negative. Also the second moment violates 
scale invariance and is negative if b < 1 . 

Given a system of points Xj we can approximate the singular integral 
by a sum 



^ /''/('')$ £<<•,/(•'•,) (C6) 



rb . , ,dx 

i=i 

where the weights are determined as above by solving the system 



= w i xk j for = 0, • • ■ z/ - 1. (C.7) 

The choice of equally spaced points gives good answers in many cases. 

These methods are used in the text (towards the end of the second chap- 
ter) to solve the integrals equations for the wavefunction of the baryon. 
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